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Abstract. Changes in the hydrologic balance in many irri- 
gation areas, including those in the Murray Basin, Aus- 
tralia, have resulted in high watertables and salinity prob- 
lems. However, where suitable aquifers exist, ground- 
water pumping and subsequent irrigation application af- 
ter mixing with surface waters (referred to as conjunctive 
water use) can control salinity and watertable depth and 
improve productivity of degraded land. In order to assess 
where conjunctive water use will successfully control 
salinity, it is necessary to estimate the effects of pumped 
groundwater salinity on rootzone salinity. A simple 
steady rate model is derived for this purpose from mass 
conservation of salt and water. The model enables an 
estimate to be made of rootzone salinity for any particu- 
lar salinity level of the groundwater being used in con- 
junction with surface water; this enables calculation of 
the required crop salt tolerance to prevent yield reduc- 
tions. The most important input parameters for the 
model are groundwater salinity, the annual depth of class 
‘A’ pan evaporation, the annual depth of rainfall, the 
salinity of irrigation water, and a leaching parameter. For 
model parameters nominated in this paper, where 
groundwater salinity reaches 5 dS/m a crop threshold salt 
tolerance greater than 1.6 dS/m is required to avoid yield 
reductions. Where groundwater salinity approaches 
10 dS/m, a crop threshold tolerance of 3 dS/m is requir- 
ed. Whilst the model derived indicates that rootzone 
salinity is sensitive to groundwater salinity, rootzone 
salinity is insensitive to leaching for leaching fractions 
commonly encountered (0.1 to 0.4). The insensitivity to 
leaching means that it could be expected that similar 
yields could be attained on heavy or light textured soils. 
This insensitivity also implies that there is no yield penal- 
ty from increasing the mass of pumped salt by pumping 
to achieve maximum watertable control in addition to 
leaching. The model developed is also used to estimate 
yield reductions expected under conjunctive use, for any 
particular levels of groundwater salinity and crop salt 
tolerance. 


Correspondence to: B. Prendergast 


In inland semi-arid areas, irrigation can result in large 
increases in agricultural productivity. However, if natural 
groundwater flows are restricted, high watertables and 
insufficient leaching can prevail. Where unconfined 
aquifers exist, groundwater pumps can provide effective 
drainage. Sub-surface drainage problems are often exac- 
erbated if there are no nearby rivers or oceans to provide 
economic export of saline effluent from groundwater 
pumps or tile drainage. Although evaporation basins can 
be used to concentrate drainage effluent, this infrastruc- 
ture is costly and often produces adverse environmental 
conditions. Where the salinity of effluent is low there is 
potential for its irrigation application, after being mixed 
with surface waters. Conjunctive use of drainage effluent 
for irrigation has many benefits including: 

(a) obviating the immediate need for export of drain- 
age effluent, whilst maintaining effective leaching and 
therefore agricultural productivity, 

(b) reducing salt export requirements to a level neces- 
sary only to prevent groundwater degradation, 

(c) providing additional irrigation water where ground- 
water salinity is low and when accelerated groundwater 
degradation may be acceptable, 

(d) better control of rootzone waterlogging without the 
need for export of groundwater, 

(e) having diffuse dispersion of saline effluent on the 
land surface with reduced environmental problems asso- 
ciated with evaporation basins; problems which can in- 
clude the hazard of accumulation of toxic elements and 
unacceptably high levels of soil and groundwater salinity 
under and near basins. 

Conjunctive use, however, also has an impact on root- 
zone salinity due to applied salt loads. The most important 
factors affecting rootzone salinity and productivity are: 

(a) the salinity of groundwater, 

(b) the amount of groundwater applied, 

(c) the total amount of applied water, 

(d) the amount of water lost below the rootzone in- 

cluding the leaching fraction, and 

(e) the salt tolerance of the crop. 
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Toxic elements and pollutants can also severely restrict 
the potential for conjunctive use, although this aspect is 
not considered here. 

Prendergast and Heuperman (1988) suggested that the 
effect of conjunctive water use on rootzone salinity could 
be modelled through use of a leaching equation with the 
steady state equations of the salt and water budget. 
Slavich (1992) used this approach to assess the potential 
for conjunctive use; however, the groundwater budget 
equation was not considered, and the analysis was re- 
stricted to the salt and water fluxes through the rootzone; 
this meant that the model could not be applied to field 
conditions, because under field conditions of ground- 
water pumping, recharge in addition to leaching fluxes 
enters the groundwater. In addition, Slavich (1992) uses 
the somewhat arbitrary criteria of achievement of zero 
“net recharge” as an indicator of the effectiveness of con- 
junctive use for controlling waterlogging and salinity, yet 
it is not clear what the relationship is (if there is one) 
between net recharge and crop yield. 

In this paper the importance to rootzone salinity of the 
factors (a) to (e) listed above will be assessed using the 
equations of mass conservation of salt and water in a 
simple model of a conjunctive use system. The important 
relationship between groundwater salinity and rootzone 
salinity is developed as suggested by Prendergast and 
Heuperman (1988), but extended to include an uncon- 
fined aquifer, so that crop salt tolerance required to pre- 
vent yield reductions under conjunctive use can be calcu- 
lated from the model. The model is then used to assess 
management practices which can increase the effective- 
ness of conjunctive use in controlling salinity. Additional 
equations are developed to estimate reduced crop water 
use under saline conditions, so that crop yield under con- 
junctive use can also be estimated as a function of 
groundwater salinity and crop salt tolerance. In a forth- 
coming paper we will apply the model using experimental 
data from the Shepparton region. 

Steady rate equations have been shown to provide 
good prediction of water salinity effects on crop yield in 
lysimeters (Letey et al. 1985), and under field conditions 
(Prendergast 1993). In this paper the steady rate equa- 
tions used by Prendergast (1993) are extended to include 
groundwater mass balance from an unconfined aquifer. 
The only major assumption additional to those in models 
developed by Letey et al. (1985) and Prendergast (1993) is 
that of mass conservation of groundwater, because water 
salinity in the conjunctive use model is determined by the 
volume and salinity of pumped groundwater. This addi- 
tional assumption of mass conservation is used in almost 
all models without the requirement for field validation. 
Use of the validated steady rate equations means that the 
most rigorous testing for the conjunctive use model has 
already been undertaken. Experimental data to test the 
model at the scale of the area of influence of a ground- 
water pump would encompass excessive variability be- 
cause of extraneous factors that are more likely to affect 
large scale experiments, and because of practical aspects 
of sampling (including cost) that limit the amount of ex- 
perimental data that can be obtained. This greater vari- 
ability associated with large experiments means that if the 


model predictions are tested at catchment scale they are 
much more likely to fall within the range of experimental 
error. Therefore, because the steady rate equations have 
already been validated, application of the conjunctive use 
model is possible once the model parameters have been 
determined. 


Abbreviations 


A Salinity threshold (maximum value of rootzone salinity, EC,, at 
which no yield reduction is evident) in Eq. (9) of Maas and Hoffman 
1977 (dS/m); B percent yield decrease per unit rootzone salinity 
increase in Eq. (9) of Maas and Hoffman 1977 (%/dS/m); Bo yield 
response factor in Eq. (11) of Stewart et al. 1977; C, average ground- 
water salinity (dS/m); C, salinity of rainfall (dS/m); C, average 
concentration of salts in the soil solution in the rootzone (dS/m); 
C,, concentration of salts in the soil saturation extract in the root- 
zone (dS/m) (= EC,); C; average concentration of salts in total water 
applied to the crop (dS/m); C,, salinity of surface supply water for 
irrigation (dS/m); D annual depth of exported drainage ground- 
water (m/y); EC, concentration of salts in the soil saturation extract 
in the rootzone (dS/m) (= C,,); E, annual “Class A” pan evapora- 
tion (m/y); ET annual crop evapotranspiration (m/y); ET,, annual 
maximum annual crop evapotranspiration (m/y); F annual equiva- 
lent ponded depth of soil water bypass flows (m/y); G annual depth 
of upward groundwater leakage into the conjunctive use system 
(m/y); J annual total depth of water applied to the crop as irrigation 
and rainfall (m/y); J fitted coefficient in equation of Rhoades 1974 
(dimensionless); K non-leaching recharge — annual equivalent pond- 
ed depth of groundwater recharge which does not contribute to 
rootzone leaching in the conjunctive use area (equal to S+ F+G) 
(m/y); LF the leaching fraction; the fraction of applied water which 
passes below the crop rootzone (dimensionless); P annual equiva- 
lent ponded depth of pumped groundwater (m/y); R annual depth 
of rainfall infiltrated (m/y); S annual equivalent ponded depth of 
groundwater recharge from irrigation supply channels or pipe deliv- 
ery systems (m/y); W annual infiltrated depth of surface water sup- 
plied for irrigation (m/y); Y yield under salinity induced water stress 
relative to maximum yield (dimensionless). 


Theoretical development 


A conjunctive use system in an unconfined aquifer can be 
conceptualised as in Fig. 1. The model approach assumes 
steady rate processes. This assumption is well satisfied if 
sufficient time is allowed for a dynamic equilibrium to be 
achieved. It is also assumed initially that pumped ground- 
water is spread over the whole area of influence of the 
groundwater pump, because spreading groundwater over 
a large area reduces the salinity of applied water under 
conjunctive use, and therefore represents good salinity 
management. The case when this practice is not followed 
is examined in a later section of the paper. The model 
takes no account of the dissolution or precipitation of 
salts in the profile, or the extraction or excretion of salts 
by plants, these processes being considered quantitatively 
negligable compared to other terms. The assessment of 
the effect of groundwater salinity on yield is made by 
assuming groundwater salinity can take any constant 
value. This assumption is appropriate because the rate of 
groundwater degradation under conjunctive use is low, 
although the rate is influenced by a number of factors 
(Prendergast et al. 1994). The level of groundwater salin- 
ity at some future time can be estimated as outlined by 
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Fig. 1. Conceptual model of irrigation with conjunctive water use. 
Recharge to the unconfined aquifer occurs from leaching fluxes 


under crops (/- LF), channel seepage (S), upward leakage (G), and 
soil-water bypass flux (F) 


Prendergast et al. (1994). Thus the model presented here 
can also be used to predict the effects of conjunctive use 
on crop yield at future times, by using in the model the 
predicted level of groundwater salinity at any time of 
interest. 

In the analysis presented here conventional units of 
electrical conductivity are used as a surrogate for salinity 
concentration. This implies that the relationship between 
solute concentration and electrical conductivity is as- 
sumed constant for all equations, because mass conserva- 
tion is expressed in terms of conductivity. 

In Fig. 1, J is the total depth of applied water, includ- 
ing the depth of irrigation water (W), rainfall (R), and 
pumped groundwater (P). The volume of groundwater in 
the system is assumed seasonally constant and is re- 
charged from upward leakage (G), leaching fluxes under 
irrigated crops (I: LF), channel seepage (S), and soil- 
water bypass fluxes (F), under irrigated crops. The leach- 
ing fraction (LF), is defined as the fraction of applied 
water passing below the rootzone. At very low pumping 
rates the resulting restricted drainage to the aquifer 
means the leaching fraction will be principally deter- 
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mined by the pumping rate. At high pumping rates (and 
resulting free drainage to the aquifer) the leaching fraction 
will be determined principally by soil type and stratigra- 
phy, though agronomic, water management and climatic 
factors will also have an influence. 

In this paper, based on the assumption of steady rate 
processes, all rates are expressed as their annual average 
values. 

For analysis of the conceptual model illustrated in 
Fig. 1, groundwater recharge has been divided into two 
fractions which have different impacts on rootzone salin- 
ity and therefore productivity. These two fractions are: 

(a) that portion (I - LF) which leaches salts from the 
crop rootzone in the area of conjunctive use, and 

(b) recharge (K) that does not provide leaching in the 
conjunctive use area, represented by K = F+S+G (see 
Fig. 1). 


The implications of groundwater salinity for crop suitability 
in an irrigation region 


At the steady state achieved under steady rate fluxes, 
mass conservation of groundwater for the model of Fig. 1 
requires that the equivalent ponded depth of pumped 
groundwater, P, is equal to the inputs to the groundwater 
system (J - LF + K), i.e.: 


P=I-LF+K. (1) 


At equilibrium the depth of infiltrated water (/) is the 
sum of the depths of infiltrated channel water (W), infil- 
trated rainfall (R), and pumped groundwater (P), so that: 


1=R+W+P. (2) 


The equivalent statement for mass conservation of 
salt is: 


IC, = RC, + WC,+ WC, + PC,, (3) 


where C;, is the average concentration of salts in applied 
water, C, is the average concentration of salts in rainfall, 
C,, is the salt concentration of irrigation channel supply 
water, and C, is the groundwater salinity in the aquifer. 

Substituting P from Eq. (1) into (3) gives the mass of 
salt applied at the soil surface to crops, as: 


IC, = RC,+ WC, +C,(I- LF +K). (4) 


Equation (4) indicates that for any particular ground- 
water salinity (C,), other values being constant, the mass 
of salt applied to the crop will depend on the depth of 
groundwater recharge through leaching. This equation 
for the conjunctive use system can be used to provide 
management guidelines where the effect of irrigation wa- 
ter salinity on crop yield is known under local conditions. 
However, it is believed that crop yield is more functional- 
ly related to rootzone salinity than to irrigation water 
salinity, and crop salt tolerance guidelines are expressed 
in terms of rootzone salinity (see Maas and Hoffman 
1977). Thus a relationship between irrigation water salin- 
ity and rootzone salinity is now sought. 
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Whilst it is generally accepted that crop yield is affect- 
ed by both the rate of salt application and the rate of salt 
leaching, there is conflicting evidence as to the relative 
importance of these two factors. One view is that yield is 
most closely related to the mean soil salinity in the root- 
zone (Bower et al. 1969; Shalhevet and Bernstein 1968; 
Bower et al. 1970; Ingvalson et al. 1976; Rhoades 1974; 
Ayres and Westcot 1985; Hoffman et al. 1979; Prender- 
gast 1993). Another view is that crop response is more 
highly correlated with irrigation water salinity, and may 
be insensitive to high salinity in parts of the soil profile 
(Bingham and Garber 1970; Lunin and Gallatin 1965); 
this is implied by use of an equation for the water uptake 
weighted mean salinity (Dirksen 1985), derived by Bern- 
stein and Francois (1973) and suggested by Smith and 
Hancock (1986) to be an appropriate yield indicator. 
Other approaches have also been used as in Ayers and 
Westcot (1985) and Hoffman and van Genutchen (1983). 
Prendergast (1993) illustrated that both the equations of 
Rhoades (1974) and Hoffman and van Genutchen (1983) 
provided good description of field data. 

The semi-empirical leaching equation of Rhoades 
(1974) will be used for the development of the conjunctive 
use model, because of its simplicity and its close relation- 
ship with mean soil salinity, which appears to be the best 
yield indicator in many of the studies listed above. The 
equation is: 


C,=0.5J-C,[1+1/LF], (5) 


where J is an empirically determined coefficient (and in 
Eq. (5) replaces K in the original equation of Rhoades 
(1974) to avoid symbol duplication). Rhoades (1974) sug- 
gests J = 0.8 approximately. In Eq. (5) C, can be related 
to rootzone salinity measured on the soil saturation ex- 
tract (C,.), which is used in salt tolerance guidelines, 
through C, = 2C,., by assuming the water content of the 
saturation extract is about twice that at field capacity 
(Maas and Hoffman 1977; van Schilfgaarde et al. 1974). 

In the development of the conjunctive use model the 
depth of applied water in the mass balance equation (4) 
can be replaced by a climatic factor to produce a more 
general expression. Because the volume of groundwater 
in the system of Fig. 1 is assumed constant, yearly crop 
evapotranspiration, ET , is the sum of total applied water 
depth, J, minus the water lost to the groundwater 
(1: LF + F), that is: 


ET=I(1-LF)—F. (6) 


Soil-water bypass flows in the Shepparton region in the 
Murray Basin, Australia, are estimated to represent up to 
8 mm/yr on a major soil type, which is negligibly small 
compared with other terms in Eq. (6). Equation (6) will 
therefore be simplified to ET = I(1— LF). It may not be 
possible to ignore F in the groundwater budget equation 
(1), because the magnitude of water fluxes is much smaller 
than in Eq. (6), and F can be a substantial component in 
the groundwater budget. 

When maximum productivity is achieved, crop evapo- 
transpiration approaches maximum evapotranspiration 
(ET,,) (Doorenbos and Kassam 1973). Maximum crop 
evapotranspiration, and therefore maximum yield, is usu- 


ally aimed for in irrigation areas because of the high cost 
of irrigation layout. ET,, can be approximated by the 
evaporation rate (E,) from a well sited evaporation pan 
with a pan coefficient of 0.85, and a crop factor of 1 
(Doorenbos and Pruitt 1977) as: 


ET,, = 0.85E, . (7) 


Taking ET = ET,, and F = 0 in Eq. (6), with J =0.8 
and C, =2C,, in Eq. (5), Eqs. (1), (2), (3), (5), (6) and (7) 
yield: 


1 
C,, =0.2 (: + iF) (8) 


K(C,~C 
iy 0.85E, 


(1—LF)+LF(C,— ca 


Equation (8) describes the effect of conjunctive use of 
groundwater on rootzone salinity for the system in Fig. 1 
for full irrigation, LF being dependent on P (Eq. (1)). 
C,, is plotted against the leaching fraction in Fig. 2 for 
K=0.06 m/yr. 

Figure 2 illustrates the mean crop rootzone salinity 
expected under irrigation incorporating conjunctive water 
use for certain nominated values of the model parameters. 
As shown in Fig. 2, C,, is very sensitive to restricted leach- 
ing (LF < 0.05) but is insensitive to leaching fractions 
commonly encountered (LF from 0.1 to 0.4). Light soils 
lead to higher LF and therefore higher P. The beneficial 
effects of higher LF under conjunctive use are offset by 
the greater depth of pumped groundwater which leads to 
higher C,;. Rootzone salinity increases only slightly for 
very high LF (e.g. LF > 0.3) because of the greater depth 
of pumped groundwater (and therefore greater mass of 
surface applied salt). Figure 2 also indicates that C,, is 
more sensitive to LF at higher values of C,. 


C,=20 dS/m 


C,=10 dS/m. 


dS/m 
= 2 dS/m 


Rootzone Salinity, C,, (dS/m) 


T 


0.1 0. 2 0.3 0. 4 0.5 


Leaching Fraction, LF 


Fig. 2. Rootzone salinity (C,,) under conjunctive water use for max- 
imum crop evapotranspiration. Plot of Eq. (8) with C, = 0.1 dS/m, 
K = 0.06 m/y, R = 0.45 m/y, C, = 0.008 dS/m, E, = 1.4 m/y 
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The highest value of rootzone salinity at which no 
yield reduction occurs is defined as crop threshold salinity 
(A). For a threshold value of 1.6 dSm~! full yields would 
be expected to be maintained for C, below 5dSm“', at 
least for LF between 0.15 and 0.5 (see Fig.2). A crop 
threshold salinity of 3.0 dS m~' would be required for C, 
up to 10dSm“', provided LF > 0.1. 


The effect of groundwater salinity on crop yield 


In many irrigation regions with salinity problems it is 
likely that some crops will suffer yield reductions and 
lower water use, but can still be grown economically. In 
order to assess management strategies under such condi- 
tions, the model must include crop water use estimated 
for rootzone salinities greater than crop threshold values. 

Maas and Hoffman (1977), suggest that relative yield, 
Y (percent of maximum yield, Y,,,,), can be expressed as 
a function of average rootzone salinity, that is: 


% Yield = 100 — B(EC, — A), EC,>A, (9) 


where A is the salinity threshold, B the percent yield de- 
crease per unit salinity increase and EC, is the average 
salinity of the soil saturation extract in the rootzone. In 
dimensionless form, after substituting C,, for EC, to elim- 
inate differences in notation, and expressing yield in frac- 
tional form (i.e. Y’= Yield/Y,,,,), Eq. (9) becomes: 


Y’=1—0.01B(C,,— A), C,,>A. (10) 


Decreases in osmotic potential of soil water, associat- 
ed with high rootzone salinity, reduces crop water use 
(Bresler et al. 1982; Bernstein 1975), which is highly corre- 
lated with crop yield (de Wit 1958; Rasmussen and Hanks 
1978; Shalhevet et al. 1976). Crop yield can be expressed 
as a function of crop evapotranspiration and maximum 
evapotranspiration (Doorenbos and Kassam 1979; Stew- 
art et al. 1977; Hanks 1983) by: 


(11) 


where f, is a yield response factor (Doorenbos and Kas- 
sam 1979; Hanks 1983). Solution of Eqs. (7) and (11) 
yields the expression for evapotranspiration by crops ex- 
periencing salinity induced water stress: 


ET =0.85E,(Y¥ —1 + Bo)/Bo - (12) 


An expression for relative crop yield (Y) can then be 
found using Eq. (12) instead of (7) for crop water use, and 
solving Eqs. (1), (2), (3), (5), (6), (10) and (12), with F = 0 in 
Eq. (6) and J = 0.8 in Eq. (5), to yield: 


0.2 K(C,— Cy) C,) 
Y=1—0.01B {(02+ 75) 


“(1 


Equation (13) has been illustrated in its simplest form 
above, by leaving relative yield (Y) on the right hand side 
as well as the left. The equation is only valid for rootzone 
salinity (C,.) above the threshold (A) (ie. Y < 1). 


1.0 7 
C,= 5 dS/m 
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Fig. 3. Yield of fully irrigated perennial pasture suffering salinity 
impaired water use. Plot of Eq. (13), C,. > A, with C, = 0.1 dS/m, 
R=045m/y, C,=0.008dS/m, K=0.06m/y, E,=14m/y, 
A= 1.6 dS/m, B= 8.9 and B, =1 


The quantitative implications of this theory can be 
illustrated by substituting values for parameters. These 
values are chosen as representative for perennial pasture 
grown in the Shepparton region in the Murray Basin, 
Australia, where channel supply water is generally of high 
quality, C, = 0.1 dS/m being typical. Using an average 
regional rainfall R = 0.45 m/yr, C, = 0.008 dS/m, E, = 
1.4 m/yr and K of 0.06 m/yr, Y is plotted as a function of 
the leaching fraction and the groundwater salinity in 
Fig. 3. In this example crop salt tolerance parameters are 
those measured by Mehanni and Repsys (1986) for peren- 
nial pasture, giving A = 1.6 and B = 8.9 in Eq. (9). Equa- 
tion (13) is valid for 12.9 > C,, > A. The upper limit for 
C,, is represented by zero yield in Eq. (9). In reality Eq. (9) 
loses linearity at very high values of rootzone salinity, and 
these levels are economically unrealistic. A value of 1.0 is 
assumed for the yield response factor (8), although this 
can be substituted for a more appropriate value, where 
known. The yield response factor (f,) often has a value of 
about 1 (Doorenbos and Kassam 1979) and has a value of 
about 1.05 for lucerne (Hanks 1983). Equation (13) has 
been solved iteratively and plotted in Fig. 3 using these 
values. 

From Fig. 3 it can be seen that for LF > 0.2 approxi- 
mately when C, = 10 dS/m, an increase in the leaching 
fraction results in increased yield. This higher yield with 
greater leaching becomes more important in absolute 
magnitude with higher levels of groundwater salinity. At 
low levels of groundwater salinity, yield is not particular- 
ly sensitive to the range of leaching fractions typical of 
many soils (i.e. 0.1—0.4). 

An estimate of crop yields under conjunctive use, for a 
particular level of groundwater salinity, can be made 
from Fig. 3. For example where groundwater salinity is 
10 dS/m, conjunctive water use will result in 88% of full 
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yield. Similarly a groundwater salinity of 8 dS/m will re- 
sult in 92% of maximum yield. 

Reduction in leaching to low values of LF causes de- 
creased yield for any particular level of groundwater 
salinity (Fig. 3). As leaching approaches zero, rootzone 
salinity increases indefinitely (Fig. 2). In Fig. 3, however, 
the asymptote of the equation shifts to higher values of 
LF with increases in C,, so that exceeding a minimal 
leaching becomes more critical in comparison with Fig. 2. 
This occurs because K, assumed constant, represents a 
greater percentage of the water budget when the crop uses 
less water whilst under salinity stress. 


Salt export for salt balance 


Under irrigation with conjunctive water use, ground- 
water is cycled, so that salt imports in irrigation water 
and rainfall can cause groundwater degradation. Salt ex- 
port can maintain a regional salt balance if this export 
balances surface salt import (Prendergast et al. 1994). The 
depth of groundwater used conjunctively is then reduced 
by export of groundwater out of the system, and rootzone 
salinity is lower. Export in this paper refers to discharge 
of drainage effluent from the irrigated area through water 
conveyance infrastructure on the surface. 

Mass conservation of salt for the system with salt ex- 
port means that the annual applied salt (JC;), is made up 
of salts in rainfall (RC,), plus salts in irrigation supply 
water (WC,,), and pumped groundwater (PC,), minus the 
depth of exported groundwater (D) times its salinity (C,). 
D must be less than the depth of groundwater recharge, 
I-LF+K (or P, through Eq. (1)), otherwise a falling 
watertable results and steady state assumptions used for 
the derivation of the conjunctive use model are not valid. 
The equation is therefore: 


IC, = RC,+WC,+(P—D)C,, P>D. (14) 


Salt balance maintained through export of some 
groundwater means that the mass of salt export, DC,, 
equals the salt mass applied in rainfall and irrigation 
water, 1.e.: 


DC,=RC,+WC,. (15) 


If D = P, then by definition there is no conjunctive use. 
From Eq. (15) the depth of groundwater for export (D) is 
(RC, + WC,)/C, and this equals P (or LF + K from 
Eq. (1)), when C, =(RC, + WC,)/(I- LF +K). Thus the 
steady state condition that P > D requires that: 


C, > (RC, + WC,)/(I- LF + K). (16) 


The relationship between rootzone salinity, irrigation 
water salinity and the leaching fraction in Eq. (5) in the 
presence of groundwater export and with J = 0.8 can be 
solved with (6), (7), (14), (15), and (16), with F = 0 in (6), to 
give: 


1 
=0.2C, (1+ — 
+ (uF + 


K(1—LF) 
O.85E, 


RC, + WC, 
I-LF+K 


(17) 
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Fig. 4. Rootzone salinity under conjunctive water use for maximum 
crop evapotranspiration and groundwater export for salt balance. 
Plot of Eq. (17) with E, = 1.4 m/y and K = 0.06 m/y 


Equation (17) is plotted in Fig. 4 using the same pa- 
rameter values typical of the Shepparton Irrigation Re- 
gion in Australia as were used in Fig. 2. When compared 
with Fig. 2, Fig. 4 (which also describes crops under no 
salinity stress) indicates that groundwater export to 
maintain salt balance, described by Eq. (15), has little ef- 
fect on rootzone salinity. This is because D (as determined 
from Eq. (15)) is only a small fraction of P, except at very 
low levels of groundwater salinity, in which case rootzone 
salinity is well below threshold in the conjunctive use 
system. 


Discussion 


Management to minimise rootzone salinity 
under conjunctive use 


The effect that groundwater recharge from the surface has 
on the plant rootzone in Fig. 1 is dependent on the path 
these water fluxes take to the watertable. Water move- 
ment through the soil matrix has two opposing effects on 
rootzone salinity in the conjunctive use system: 

(a) this water leaches salts from the crop’s rootzone, 
thus lowering rootzone salinity, but 

(b) these leaching fluxes increase groundwater recharge 
and hence the depth of groundwater that has to be 
pumped at equilibrium, and therefore indirectly increase 
rootzone salinity. 

Non-leaching recharge (K) affects the plant rootzone 
only by increasing the depth of pumped groundwater. 

Water movement through the soil matrix is generally 
low when heavy textured soils are irrigated and high 
when light soils are irrigated. Water application depths, 
irrigation frequency, and the rate of groundwater pump- 
ing also affect water movement through the soil matrix. 
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However, the above analysis has shown that for the 
parameter values adopted (at least), changes to water 
movement below the rootzone through the soil matrix, as 
indicated by LF, have relatively little impact on rootzone 
salinity for LF > 0.1. Pumps should therefore be operat- 
ed to achieve free drainage to the aquifer in order to 
ensure sufficient leaching. Pumping can then also lower 
the watertable, which will improve crop yields where 
rootzone waterlogging is a significant problem. Increas- 
ing the pumping rate, however, increases well drawdown, 
which also has the potential to increase G and S$ and 
therefore K. If rootzone waterlogging is not a significant 
problem, pumping at lower rates is therefore desirable to 
reduce pumping costs. Pump management guidelines, 
therefore, must be evaluated in the context of local condi- 
tions. 

The importance of reducing non-leaching recharge for 
different given aquifer salinities can be assessed by plot- 
ting rootzone salinity against K in Eq. (8). Figure 5 pre- 
sents the results of these calculations. It is notable in 
Fig. 5 that under the ideal conditions represented by 
K =0, groundwater of much higher salinity can be 
utilised without excessively increasing rootzone salinity. 
For example, with a groundwater salinity of 10 dS/m, 
LF =0.1 and a crop threshold salt tolerance of 2.3 dS/m 
there would be little yield reduction (Fig. 5). 

Whilst rootzone salinity is relatively insensitive to the 
leaching fraction over the higher range of LF commonly 
achieved (as shown in Figs. 2 and 4), C,, is linearly sensi- 
tive to non-leaching recharge. Thus considerable yield 
benefits can be obtained from reduction in non-leaching 
recharge. For example, it can be seen from Fig. 5 that a 
reduction in non-leaching recharge from 0.25 m/yr to 
0.1 m/yr will result in a 35% reduction in rootzone salin- 
ity for a groundwater salinity of 5dS m7‘ and LF = 0.1. 
The benefit of reducing K is greater at higher levels of 
groundwater salinity, and greater for small values of LF 
(Fig. 5). Conversely, the detrimental affect of greater ap- 
plied salt loads from increased K can be partially offset by 
greater rootzone leaching. The minimum value of root- 
zone salinity as a function of LF in the solutions to Eq. (8) 
plotted in Fig. 2, is also implied in Fig. 5. It can be seen 
from Fig. 5 that for low values of K, increasing leaching 
from LF = 0.01 to LF = 0.4 actually increases rootzone 
salinity. At high values of K this addional leaching re- 
duces rootzone salinity. Figure 5 shows that the value of 
K at which higher leaching begins to reduce rootzone 
salinity is greater with higher levels of C,. 

Management of the conjunctive use system can reduce 
K in several ways. Reduced channel seepage (S) can be 
achieved in some regions by piping irrigation water. A 
reduction in channel supply water salinity can reduce 
seepage by reducing hydraulic conductivity. Water appli- 
cation with sprinklers is may offer a reduction in bypass 
flows (F). 

If the depth of pumped groundwater used for irriga- 
tion exceeds groundwater recharge in a conjunctive use 
system, the extra irrigation groundwater must have been 
made up of an equivalent mass of groundwater flowing 
laterally into the system from outside, if the system is in 
steady state. Such lateral flows can be considered as 
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Fig. 5. Rootzone salinity vs. non-leaching recharge for full irrigation 
with conjunctive water use with no salt export. Plot of Eq. (8) with 
C,, = 0.1 dS/m, R = 0.45 m/y, C, = 0.008 dS/m, E, = 1.4 m/y 


adding to K, since they do not provide leaching in the 
conjunctive use system. Although K, as defined earlier, 
does not describe the situation where groundwater used 
conjunctively is irrigated over only part of the pump’s 
area of influence, the term K could be readily modified so 
that this would be included. This situation arises if the 
area to which conjunctively used groundwater is applied 
is smaller than the area of influence of the groundwater 
pump, thus providing a greater depth of groundwater for 
crops irrigated on that reduced area. However, in order to 
minimise non-leaching recharge, it is important that con- 
junctively used groundwater be irrigated over a large per- 
centage of the pump’s area of influence. In a conjunctive 
use system the depth of groundwater recharge is small in 
comparison with the depth of crop water use. Therefore 
use of groundwater to provide a large proportion of crop 
water use implies that much of this groundwater has 
come from outside the area of conjunctive use. Minimis- 
ing K and thus rootzone salinity under conjunctive use 
therefore precludes a large proportion of crop water use 
being provided by pumped groundwater. 

The effectiveness of the conjunctive use system for con- 
trolling salinity also depends on management of the 
groundwater catchment, of which the conjunctive use sys- 
tem is a part. Regional hydrogeological characteristics 
may result in high values for G (Fetter 1988), so that 
groundwater pumping in a conjunctive use system may 
result in beneficial leaching of upstream areas where there 
are no groundwater pumps. A policy of conjunctive water 
use can allow crop threshold salt tolerance properties to 
be exploited to enable control of irrigation salinity with- 
out loss of crop yield. The effectiveness of conjunctive use 
over a catchment is therefore maximised by utilising 
threshold tolerance over as large an area of the catchment 
as possible. Irrigated areas which have leaching but not 


173 

10 

8 

4 

f 

0 
1 

| 


174 


conjunctive use, and hence little salt application, do not 
utilise crop threshold salt tolerance. When the ground- 
water catchment is a closed groundwater system, leached 
areas without conjunctive use can cause high values of K 
in conjunctive use areas which overlie groundwaters 
downstream. Surface redistribution of groundwater back 
to upstream leached areas which have no installed pumps 
increases the area over which crop threshold tolerance is 
utilised for disposal of groundwater. Thus redistribution 
can improve the effectiveness of conjunctive water use in 
controlling salinity in the catchment. Surface redistribu- 
tion of groundwater is most important in groundwater 
catchments where hydrogeological characteristics result 
in high values of G in the conjunctive system in parts of 
the catchment and low or negative values of G in other 
parts. 


Rootzone leaching 


As stated earlier, increasing LF has two opposite effects 
on crop rootzone salinity in a conjunctive use system; 
firstly, increased LF results in greater salt transport out of 
the rootzone, which reduces rootzone salinity, but se- 
condly, by resulting in the need for a greater depth of 
pumped groundwater (Eq. (1)), it increases the applied 
salt load. This means that in the conjunctive use equa- 
tions (8), (13) and (17), there is a minimum rootzone salin- 
ity which can be achieved. This minimum represents the 
optimum leaching management of the system, and an 
expression for this state can be found by evaluating LF 
where the equation derivative is zero. Thus LF = LF,,, 
when: 


0. (18) 


The solution to Eq. (18) is a function of non-leaching 
recharge, K, and groundwater salinity, C,, so that opti- 
mal management of this system is site specific. As is illus- 
trated in Figs. 2 and 4, to the left of the minimum repre- 
sented by Eq. (18), decreased leaching will result in in- 
creases in rootzone salinity. Especially at low values of 
C,, the value of C,, is close to the minimum over a range 
of values of LF. To the right of the minimum, decreasing 
leaching will have a beneficial effect on rootzone salinity. 
An example of how LF,,, increases with K and C, is 
implied in Fig. 5. Increasing leaching from LF = 0.1 to 0.4 
increases C,, when K = 0, but decreases C,, at high values 
of K; the value of K where C,, is equal for these two 
leaching fractions and a particular C,, increases with C, 
(see Fig. 5). 

From the model describing the conjunctive use system 
illustrated in Fig. 1, the value of LF,,, depends on the 
importance of leaching in controlling rootzone salinity. 
In irrigation areas where the parameters used in Figs. 2 
and 3 are valid, LF,,, occurs at approximately LF = 0.35. 
For areas where leaching is of lesser importance in con- 
trolling rootzone salinity and the equation of Bernstein 
and Francois (1983) is valid, for example, LF,,, occurs at 
low leaching fractions, and reduction in leaching will im- 
prove rootzone salinity for values of LF commonly en- 


countered. This highlights the importance of local deter- 
mination of the leaching parameters before application of 
the conjunctive use model. 

Equation (5) can be fitted to a range of different soils 
through its empirically determined coefficient, J. Root- 
zone salinity as described in Eqs. (8) and (17) is propor- 
tional to J, so that a higher value of J means greater 
rootzone salinity in the conjunctive use model. The deriv- 
ative OC,./0(LF) for Eqs. (8) and (17) is also directly pro- 
portional to J, so that the gradient of the lines in Figs. 2, 
3 and 4 increases with increasing J. 


Conclusion and comment 


Where the leaching equation (5) of Rhoades (1974) is 
valid, the conjunctive use model developed here suggests 
that rootzone salinity is insensitive to leaching in the 
range of leaching fractions commonly encountered (from 
0.1 to 0.4). Strategies such as intermittent pumping, which 
can reduce the pumped water depth and thereby min- 
imise leaching, could result in decreased crop yield if low 
leaching results (e.g. LF > 0.05). The conjunctive use 
model indicates that high induced leaching rates will be 
unlikely to increase rootzone salinity. In areas where 
rootzone waterlogging is a problem, lowering the water- 
table could therefore be achieved by using high pump- 
ing rates or more continuous pumping to increase the 
pumped water depth, if this does not increase K exces- 
sively. 

The model indicates that where the nominated param- 
eters are applicable, and the groundwater salinity is 
10 dS/m, a threshold crop salt tolerance of 3 dS/m is re- 
quired to prevent yield reduction under conjunctive 
water use. Where groundwater salinity is 5 dS/m, a crop 
threshold salt tolerance of 1.6 dS/m is required to avoid 
reductions in yield. Export of a mass of salt equal to that 
applied in irrigation water and rainfall has little effect on 
rootzone salinity under conjunctive use. 

In the analysis of the conjunctive use system, non- 
leaching recharge is found to be the most important com- 
ponent of the water budget to affect rootzone salinity. 
There is potential for large additions to this component 
when groundwater is not irrigated over the whole of the 
area of influence of the groundwater pump. In some irri- 
gation areas regional hydrogeological characteristics can 
result in excessive upward leakage (G) into the pumped 
aquifer and cause large additions to non-leaching re- 
charge. Surface redistribution of pumped groundwater to 
upstream irrigation areas without groundwater pumps 
can help counteract the adverse salinity impacts of ex- 
cessive upward leakage (G) by reducing the volume of 
groundwater to be applied in the pumped area. 

The lumped parameter model developed in this paper 
has been used to assess management options under con- 
junctive water use. The model is also likely to have appli- 
cation for providing salinity management guidelines, be- 
cause of the difficulty, and therefore cost, in obtaining 
sufficient data for an equivalent distributed parameter 
model. Applicability of the model will be enhanced by 
determination of suitable parameters through site investi- 
gations in an irrigation region. 
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Abstract. The long term sustainability of conjunctive wa- 
ter use for controlling irrigation salinity is affected by 
increase in groundwater salinity over time. This paper 
uses mass conservation of salt and water to assess 
groundwater degradation over long time scales. Manage- 
ment options which affect this rate of degradation are 
also examined. The groundwater model developed is il- 
lustrated using data from the Shepparton Irrigation Re- 
gion in the Murray Basin, Australia. The model predicts 
rapid groundwater deterioration when conjunctive use is 
conducted over only a fraction of the area of influence of 
a groundwater pump. Where the pumped aquifer is un- 
derlain by deeper groundwaters, the rate of groundwater 
degradation is also affected by leakage into or out of the 
conjunctive use system. Surface redistribution of ground- 
water from pumps installed in zones of regional ground- 
water discharge to areas recharging the regional ground- 
waters, reduces excessive degradation in the zones of 
discharge. With optimal surface distribution of ground- 
water, the rate of degradation is low. The rate of ground- 
water degradation also depends on salt inputs from irri- 
gation water and rainfall, and the average depth from the 
soil surface to the base of the aquifer. The rate of degra- 
dation resulting from applied salts in surface water and 
rainfall is typically about 0.01 dSm™~! per year for shal- 
low aquifers in the Shepparton region, but the rate is 
lower where deeper aquifers are pumped. Partial irriga- 
tion also reduces the rate of degradation because of the 
reduced rate of salt inputs. Where poorer quality ground- 
water lies within the area of influence of the groundwater 
pump, a greater rate of deterioration in the quality of 
pumped groundwater can be expected from groundwater 
mixing. In some irrigation regions limited export of 
groundwater through surface water conveyance struc- 
tures to a river is possible, so that a regional surface salt 
balance could be maintained. However, salt exports 
made equal to the rate of surface imports into the irrigat- 
ed area will only significantly impact groundwater salini- 
ty in the very long term, or where only shallow aquifers 
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can be pumped. In addition, this export can be costly for 
downstream water users, or if construction of additional 
conveyance infrastructure is extensive; export can have a 
detrimental impact on riverine ecosystems. Other man- 
agement options such as the depth of pump installation 
and the spatial distribution of irrigation water and 
pumped groundwater, which affect the redistribution of 
salts within the groundwater system, have the potential to 
have a much greater impact on local groundwater salini- 
ty. 


Where hydrogeological conditions are such that insuffi- 
cient leaching for irrigated agriculture occurs, artificial 
sub-surface drainage is necessary. If suitable aquifers ex- 
ist, groundwater pumps can provide the most economic 
means of sub-surface drainage. Application of pumped 
groundwater in conjunction with surface supplied irriga- 
tion water provides the most economic and environmen- 
tally acceptable means of dispersing drainage effluent, 
provided the groundwater is of sufficiently low salinity 
and toxicity (Salinity Program Advisory Council 1989). 
However, groundwaters underlying irrigation areas can 
have annual salt imports from irrigation water of 0.1 to 
10 tonnes per year per hectare, and this salt application 
can cause an increase in groundwater salinity through 
time. The rate of change of groundwater salinity also 
depends on whether groundwater flows result in a net 
addition or loss of salt from the conjunctive use system. 
Because conjunctive water use results in higher levels of 
salinity of applied water, sustainability of crop yields can 
also be affected by salinity and sodicity effects on soil 
physical properties, although this aspect is not addressed 
here. 

Land and water management have a large impact on 
groundwater systems, particularly where groundwater is 
pumped. For example, restricting the area of conjunctive 
use to less than the area of influence of a groundwater 
pump (referred to as spatially constricted conjunctive 
use), or reduction in the depth of irrigation water applica- 
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tion (referred to as partial irrigation), both affect the rate 
of salt accumulation under conjunctive use. Therefore the 
long term implications of conjunctive use of groundwater 
for salinity control will be dictated by groundwater man- 
agement, which controls the distribution of groundwater 
salts above and below the surface, as well as by land and 
surface water management. Good salinity management 
aims to minimise rootzone salinity and sustain productiv- 
ity and groundwater quality. The likely consequences of 
alternative management options on groundwater quality 
are examined using a simple steady rate lumped parame- 
ter model derived in this paper from mass conservation 
considerations. Lumped parameter models are usually 
used where time, but not space, is the important variable 
(Domenico 1972). Prendergast et al. (1994) used a model 
of steady rate water fluxes to determine the effect of a 
particular level of groundwater salinity on crop yield un- 
der conjunctive use; however they did not look at the 
effect of these steady rate fluxes on the change in ground- 
water salinity over long time scales, which is considered 
here. Data used in the model are taken as typical of the 
Shepparton Irrigation Region in the Murray Basin, Aus- 
tralia. 


Abbreviations 


A, Area of influence of a groundwater pump in an unconfined 
aquifer (m7); a initial groundwater salinity (Cy) minus final steady 
state groundwater salinity (c) (dS/m); b ratio of groundwater dis- 
charge out of the conjunctive use system to groundwater storage in 
the system (y~'); C, average groundwater salinity from the soil 
surface to the base of an unconfined aquifer (dS/m); C,, average 
groundwater salinity (dS/m) from the soil surface to the base of an 
unconfined aquifer within the area of influence of a groundwater 
pump, but outside the area of conjunctive water use; C,, average 
groundwater salinity from the soil surface to the base of an un- 
confined aquifer under irrigation with conjunctive water use (dS/m); 
Cy initial average water salinity (dS/m) below the soil surface; 
C, salinity of rainfall (dS/m); C,, salinity of surface supply water for 
irrigation (dS/m); c final steady state groundwater salinity (dS/m); 
E,, annual “Class A” pan evaporation (m/y); ET,, annual maximum 
crop evapotranspiration (m/y); G annual equivalent ponded depth 
of groundwater leakage (positive downwards) below the base of the 
aquifer to deeper groundwater (m/y); H depth from the soil surface 
to the base of the aquifer (m); J annual total depth of water applied 
to the crop as irrigation and rainfall (m/y); K, crop factor (dimen- 
sionless); K,, pan coefficient (dimensionless); LF the leaching frac- 
tion; the fraction of applied water passing below the crop rootzone; 
P annual equivalent ponded depth of pumped groundwater (m/y); 
R annual depth of rainfall (m/y); S annual equivalent ponded depth 
groundwater recharge from seepage from irrigation supply channels 
or pipe delivery systems (m/y); t time in years; U land fraction main- 
taining conjunctive water use within the area of influence of a 
groundwater pump; W annual infiltrated depth of surface water 
supplied for irrigation (m/y); W, annual depth of channel (or piped) 
supply water infiltrating outside the area of conjunctive water use 
but within the area of influence of the groundwater pump (m/y); W, 
annual depth of channel (or piped) supply water infiltrating within 
the area of conjunctive water use (m/y); 0 average volumetric water 
content of sediments from the soil surface to the base of the aquifer 
(m? 


Assumptions and background for theoretical development 


In deriving this model it will be assumed that groundwa- 
ter recharge into the conjunctive use system completely 
mixes with the volume of groundwater under consider- 


ation and that there is no lateral flow between the system 
and its surrounds. Since both the velocity of groundwater 
recharge and the distances along the groundwater stream- 
lines to the pump screen are highly variable, a high degree 
of mixing is expected. Even poor mixing of groundwater 
recharge with groundwater will not affect mass conserva- 
tion of salts over long time scales, although short term 
variations in groundwater salinity may occur which are 
not accounted for in this model. 

Diffusion is neglected in the model presented because 
it is assumed to be of little importance compared to the 
convective transport of solutes through leaching and 
pumping. The processes of dissolution and cation ex- 
change are also ignored as the groundwater is considered 
to be in a state of equilibrium with precipitated and ex- 
changeable ions. The rate of dissolution depends on the 
chemical nature and rate of mineral decomposition of the 
stratigraphy. Dissolution is not considered in this paper, 
although it may need to be incorporated when assessing 
irrigation regions where substantial mineral decomposi- 
tion occurs over the very short geologic time scales con- 
sidered here. 

Use of conventional units of electrical conductivity for 
salinity, in expressing mass conservation of solutes, im- 
plies that the relationship between electrical conductivity 
and solute concentration is assumed constant. 

The most important factors that cause changes in 
groundwater salinity under conjunctive use, which are 
discussed in this section, are: 

(a) salt inputs to the groundwater system from irriga- 
tion water and rainfall, 

(b) the spatial relationship between the conjunctive use 
system and the regional groundwater, and the resulting 
groundwater leakage into or out of the system, 

(c) spatial constriction of the area of conjunctive water 
use, 

(d) mixing of groundwaters, and 

(e) salt export from the irrigated area; export in this 
paper refers to groundwater being discharged from the 
irrigated area through surface conveyance structures. 

The lumped parameter groundwater model is devel- 
oped for conjunctive use in irrigation areas from the con- 
ceptual model illustrated in Fig. 1. The model is used to 
assess factors (a), (b), (c) and (e) above, but not (d), because 
the model assumptions of perfect mixing render it inap- 
propriate for analysing mixing. All water inputs and out- 
puts are based on the assumption of steady rate processes, 
with rates expressed as average annual depths. 

Crops in Fig. 1 are irrigated with surface water of 
depth W and concentration C,,, and pumped ground- 
water, P, to supplement rainfall of depth R and concen- 
tration C,. Recharge to the groundwater occurs from 
channel seepage (or pipe leakage) (S), and leaching fluxes 
below the crop rootzone (I - LF), where I is the total 
depth of applied water (equal to W + P+ R), and LF is the 
leaching fraction (the fraction of applied water passing 
below the rootzone). The magnitude of soil-water bypass 
flows which could enter the groundwater by bypassing 
the crop rootzone through large cracks and voids is as- 
sumed to be zero, which is a reasonable assumption for the 
soil type prevalent in the Shepparton Region (see Pren- 
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Fig. 1. Conceptual model of a conjunctive use system showing 
boundary conditions, inputs and outputs, for the development of a 
groundwater model for an irrigation area with conjunctive use of 
groundwater 


dergast and Noble 1990). Furthermore, bypass flows can 
be assumed to have little effect on groundwater degrada- 
tion in the conjunctive use system, because of their low 
salinity. 

Groundwater is pumped from an unconfined (or 
“watertable”) aquifer. Depth H is taken from the soil sur- 
face to the permeability layer which forms the lower 
boundary of the system. No distinction is made between 
the vadose zone and the saturated zone, so that all layers 
of the stratigraphy over H are lumped together. The aver- 
age salinity of water below the surface to depth H can 
then be C, (referred to in this paper as groundwater salin- 
ity), and the average volumetric water content over depth 
H is denoted by 0. Some leakage, G, passes through the 
low permeability layer beneath the conjunctive use sys- 
tem to deeper groundwater (Fig. 1). The area of influence 
of the groundwater pump is defined as the areal extent of 
the cone of depression of the watertable (see Todd 1980), 
as illustrated in Fig. 1. It is initially assumed that conjunc- 
tive use of groundwater and surface water for irrigation is 
maintained over the entire area of influence of the pump. 
The case where this assumption is not valid is discussed 
separately in a later section of the paper. 


H= 10m, W=0.8m 


H=100m, W=0.4m 


Groundwater Salinity, C, (dS/m) 


Time, t (Years) 


Fig. 2. Changes with time in groundwater salinity caused by salts in 
irrigation water and rainfall under conjunctive water use with G=0. 
Predictions use Eq. (2) with R=0.45 m/y, C, =0.1dS/m, C, = 
0.008 dS/m, 0 = 0.4 and C, = 4 dS/m 


Irrigation and rainfall salt inputs 


For this section the equivalent ponded depth of ground- 
water leakage beneath the system, G, is considered to be 
zero, indicating a closed groundwater system. A mass of 
salt PC, is extracted through groundwater pumping and 
applied to the soil surface for conjunctive use; this recy- 
cling of groundwater has no net effect on mass conserva- 
tion of salt within the system. Salt imports into the sys- 
tem, then, are from irrigation water (both surface applied 
and through channel seepage) and rainfall only. 

Conservation of mass of salt requires the rate of 
change of groundwater salinity (dC, /dt) to be equal to the 
annual net salt input rate divided by the volume of water 
in the conjunctive use system; i.e. 


dc, _(W+S)C,+RC, 
( 


Integrating Eq.(1) with respect to time, and using the 
initial condition C, = Cy at t = 0, gives: 
(W +S) C,+RC, 


(2) 


c= +( 


This linear relationship is plotted for two possible cases in 
Fig. 2 for parameter values typical of the range of H and 
W in the Shepparton region. An applied irrigation water 
depth of 0.8 m represents full irrigation (i.e. sufficient irri- 
gation to supply maximum crop water use) on perennial 
pastures, which utilise a large percentage of supplied irriga- 
tion water. Supply water salinity in the region is normally 
about 0.1 dS m~', although salinities up to 0.4dS/m~' 
can occur. Unconfined to semi-unconfined aquifers are 
pumped from depths from 8—150 m (Ife 1988; Nolan and 
Reid 1989; Brown 1985). More detailed description of 
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aspects of the groundwater system is given by Brown and 
Stephenson (1986); Macumber (1991); Prendergast 
(1989); Lawrence and Abele (1976) and Evans (1988). 
Figure 2 also illustrates the effect of partial irrigation, 
where water available to the crop is half the maximum 
potential crop water use in the illustrated example. Par- 
tial irrigation reduces the mass of applied salts, WC,. 
Groundwater degradation then increases with increases 
in applied irrigation water salts and with decreases in 
depth to the base of the aquifer (Fig. 2 and Eq. (2)). 


Groundwater recharge to and discharge 
from the conjunctive use system and the influence 
of the regional groundwater 


Whilst it was assumed that G=0 for Fig. 1, this is gener- 
ally not the case on average over a whole region. If indeed 
G=0 on average, then salinity problems would have been 
severe under pre-farming conditions, and it is unlikely 
that irrigation would be introduced. Despite the fact that 
regional net groundwater outflow may occur on average, 
both zones of discharge from and recharge to the regional 
groundwater from shallow aquifers will exist, depending 
on the shape of the region, boundary conditions and 
heterogeneity in the groundwater flow system (Freeze and 
Witherspoon 1967). Groundwater flows to or from the 
conjunctive use system change the mass stored and there- 
fore change groundwater salinity in the system. 

The three cases considered in this section which relate 
to groundwater flows to or from the conjunctive use sys- 
tem in Fig. 1 are: 

(1) Groundwater leakage, G, out of the system to deeper 
groundwaters. These leakage flows are usually small 
where salinity problems exist, because irrigation salinity 
problems are caused by restricted groundwater flows. In 
the Shepparton region, leakage from shallow Shepparton 
Formation aquifers to deeper aquifers has been estimated 
to range from 40 mm/y to —3.5 mm/y (Nolan and Reid 
1989). 

(2) Groundwater leakage from deeper groundwater 
upwards into the conjunctive use system (i.e. G<0). 

(3) Export of all pumped groundwater, i.e. there is no 
conjunctive use. The model developed for this purpose 
describes the effect on groundwater salinity of export to 
evaporation basins or oceans. Export may be necessary 
where groundwaters contain excessive concentrations of 
toxic elements, or where groundwater salinity is too high 
for conjunctive use. This model is also used later to de- 
scribe the situation where subsurface drainage is supplied 
by a groundwater pump, but where there is no conjunc- 
tive use on the soil surface above the drained area. 


Leakage out of the conjunctive use system 


Mass conservation of salt for the system illustrated in Fig. 1 
requires the average annual rate of change in salinity of the 
groundwater (dC,,/dt) to depend on the salt applied annu- 
ally in irrigation water and rainfall (W +S)C,+RC,), 
minus that lost through leakage to deeper groundwater 


(GC,). Thus: 


dC, _ (W+S)C,+RC, — GC, 


G>O0. (3) 
The solution to this first order linear differential equation, 
using the initial condition C,=C, at t=0, is: 


RC,+(W +S) | 

ex 
G 

RC,+(W+S)C,, 

G 


Cc. = & 


G>0. (4) 


This solution (Eq. (4)) is of the form ae~°'+c where a, b 
and c are constants. The final steady state groundwater 
salinity, c, is given by the second term on the right hand 
side of Eq. (4) (ie. c=(RC,+(W +S) C,)/G)), and so is 
determined by the ratio of the rate of salt inputs to the 
system to the rate of groundwater flow out of the system. 
Thus c is independent of the rate of groundwater recy- 
cling (ic. P). The attainment of a final steady state 
groundwater salinity implies that c must eventually be in 
equilibrium with the salinity of recharge fluxes entering 
the system, ic. the depth weighted average salinity of 
channel seepage and salinity at the base of the crop root- 
zone in Fig. 1. For a given salinity of channel supply 
water and rainfall, salinity at the base of the rootzone is 
therefore determined by G in the very long term. The 
constant a is equal to the initial groundwater salinity, Co, 
minus the final steady state groundwater salinity, c, and 
is therefore positive for groundwater salinity which de- 
creases over time, and negative for groundwater which 
degrades over time. The parameter b determines the rate 
of change of salinity to the steady state level and is given 
by the ratio of the annual depth of groundwater discharge 
out of the system to the storage in the system (i.e. G/H 6). 

Equation (4) is plotted in Fig. 3, with W set to 0.8 m/y 
to represent full irrigation for the Shepparton region. The 
depth of groundwater leakage, G, is set to 0.01 m/y, this 
value being within the range of leakage from —3.5 to 
40 mm/y, estimated for the Shepparton region (Nolan 
and Reid 1989). Other parameter values (typical of the 
Shepparton region) which are used in the evaluation of 
Eq. (4) for Fig.3 are: R=0.45 m/y, C, = 0.008 dS/m, 
C, = 0.1 dS/m, and S = 0.04 m/y. 

The exponential nature of Eq. (4) is not evident in 
Fig. 3 because of the short time scale used. For negative 
a, as illustrated in Fig. 3, the effect of leakage from the 
groundwater system is to reduce the rate of groundwater 
degradation from applied salts as described in Fig. 2, and 
establish a final steady state groundwater salinity, 
approached only after many centuries. The final steady 
state groundwater salinity, c, evaluated through Eq. (4) is 
8.6 dS/m. The rate of approach of groundwater salinity to 
the steady state value is determined by H@ and C, for 
particular values of W, C,,, R, C,, S and G. For example, 
for C, of 4dS/m, H = 30m and 6 = 0.4, Eq. (4) predicts 
that after 500 years groundwater salinity will be 5.5 dS/m 
and after 1000 years, 6.5 dS/m. The effect of these levels of 
groundwater salinity on crop yield under conjunctive use 
can be estimated as outlined by Prendergast et al. (1993). 
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Fig. 3. The impact on groundwater salinity of groundwater dis- 
charge out of the conjunctive use system. Plot uses Eq. (4) with 
R=0.45 m/y, C,=0.008 dS/m, S=0.04m, 0=0.4, H=30m, 
C, =0.1 dS/m, G=0.01 m/y and Eq. (6) after substitution of 
Eq. (8) with E, = 1.4m/y. Graphs represent full irrigation unless 
otherwise indicated 


Upward leakage into the conjunctive use system 


If G is negative (item (2) above), the rate of groundwater 
deterioration will be increased. For negative G the rate of 
salt addition to the conjunctive use system depends on the 
salinity of deeper groundwater, which is most realistically 
assumed constant over time, unlike C, in Eq. (4), so that 
Eq. (4) does not apply. For constant negative G, then, and 
constant salinity of deeper groundwater, dC,/dt is con- 
stant and C, is linear with time. As an example of this 
case, the rate of degradation can be calculated for the 
example in Fig.2 with H=10m, C,=4dS/m and 
W =0.8 m, by assuming the deeper groundwater salinity 
is equal to C,, or 4 dS/m, and the rate of leakage into the 
conjunctive use system is —0.02 m/y. This rate will add 
to the conjunctive use system a salt mass of 4dS/m x 
0.02 m x 10,000 m?/ha x 0.6 (kg/t)/(dS/m) = 480 kg/ha/y 
(using the generally accepted approximation, 1 dS/m 
= 600 mg/l = 0.6 kg/t, and assuming a water density of 
1 t/m3). This shows that a value of G= —0.02 m/y will 
double the rate of groundwater degradation in compari- 
son to the example for full irrigation given in Fig. 2. 


Salt export from the conjunctive use system 


When export to evaporation basins or to an ocean is 
available, there is usually no conjunctive water use. For 
this case deep drainage out of the system (G) is set to zero 
for simplicity, and because it is assumed that G is small in 
comparison to P. Mass conservation of salt in the con- 
junctive use system is then described by a modified 
Eq. (3), with G replaced by P, because it is now pumping 
that discharges salt from the system. For a hydrologic 
steady state, the annual depth of pumped groundwater 
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(P) must be equal to the groundwater recharge into the 
system, i.e. channel seepage, S, and the leaching flux 
under crops ((W + R) LF), so that: 


P=(W+R)LF+S. (5) 
Substituting P from Eq. (5) for G in Eq. (4) gives: 


=> & 


(W+R)LF+S 
(W+R)LF+S HO 


RC,+(W+S)C, 
(W+R)LF+S (6) 


Equation (6), an expression for groundwater salinity un- 
der conditions of complete export, is again of the form 
ae Evaluation of C, using Eq. (6) requires a 
knowledge of W. The total depth of applied water (J) 
required, consisting of crop evapotranspiration plus the 
leaching fraction can be used to evaluate W. Maximum 
crop water use can be expressed by equating crop evapo- 
transpiration with pan evaporation times a pan coeffi- 
cient (K,) of 0.85 and a crop factor (K.) of 1 (Doorenbos 
and Pruit 1977). The resulting expression for I is: 


(7) 


Since | = W+R where there is no conjunctive use, the 
depth of channel supply water (W) can be expressed as 
the difference between I, available from Eq. (7), and rain- 
fall (R). After substitution of values for K, and K, the 
expression becomes: 


0.85 E, 
(8) 


After substituting Eq. (8) for W, predictions using Eq. (6) 
are also plotted in Fig. 3 for parameter values typical of 
the Shepparton Irrigation Region. 

The effect of partial irrigation illustrated in Fig. 3 is 
obtained by replacing the crop coefficient of 1 in Eq. (8) 
with a value of 0.75, and this change results in approxi- 
mately 50% reduction in W, for R = 0.45 m. The effect of 
partial irrigation on groundwater degradation, using 
LF = 0.05 in Fig. 3, may not be simply compared with 
full irrigation for the same LF, because a reduction in 
total applied water would result in a lower value of LF. 
Because the rate of change of groundwater salinity in 
Eq. (6) depends on the coefficient of t, which is a function 
of LF and W, the lower value of LF associated with 
reduced W would tend to increase the difference between 
partial and full irrigation. 

The final steady state groundwater salinity (c) in 
Eq. (6) approaches the depth weighted groundwater 
recharge salinity into the system after a long time period. 
For c > Cy, the groundwater degrades in quality, but for 
c< Co, the groundwater improves in quality. In most 
practical cases the depth weighted recharge salinity (c) 
will be less than the groundwater salinity, as is the case for 
the evaluations of Eq. (6) illustrated in Fig. 3. If ground- 
water salinity is low and less than the depth weighted 
recharge salinity, the groundwater will probably be used 
for irrigation, in which case export is unnecessary. The 
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case where C,<c has therefore not been illustrated for 
the export case. 


Groundwater degradation from spatially constricted 
conjunctive use 


Conjunctive use of groundwater for salinity control may 
be maintained over only part of the area of influence of a 
groundwater pump. This spatially constricted conjunc- 
tive use is illustrated for an unconfined aquifer in Fig. 4. 
For simplicity the model presented here assumes water 
use by the crop is unaffected by salinity. The area of 
influence of the groundwater pump is A,, and the area 
maintaining conjunctive use UA,. The area protected 
from salinity by the groundwater pump, but without con- 
junctive use, is thus (1—U) A,. Groundwater salinity in 
the area with conjunctive use is denoted C,, and the 
depth of channel supply water applied, W,. In the area 
without conjunctive use these qualities are denoted by 
C,, and W, respectively. 

When the average salinity of recharge fluxes is lower 
than the groundwater salinity within the area of influence 
of a groundwater pump, but outside the conjunctive use 
area (Fig. 4), groundwater quality peripheral to the con- 
junctive use area in Fig. 4 will improve over time. Mass 
conservation of salt dictates that the reduction of salt 
storage in the area without conjunctive use is equal to the 
increase in salt storage in the land fraction with conjunc- 
tive use. This increase in salt storage accelerates ground- 
water degradation in the area receiving water used con- 
junctively. Groundwater degradation through spatial 
constriction of the area of conjunctive use would be ex- 
pected to become more rapid with increase in initial 
groundwater salinity, because of the greater mass of salt 
per unit of groundwater transported into the conjunctive 
use area. As the groundwater quality improves outside 
the area of conjunctive use but within the area of influence 
of the pump (i.e. the area (1 — U) A,), the rate of groundwa- 
ter degradation will slow down in the area with conjunc- 
tive use (UA)). 
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Fig. 4. Conceptual model for spatially constricted conjunctive water 
use from an unconfined aquifer 


The rate of degradation described by Eq. (9) depends on 
the applied water depth both inside and outside the area 
of conjunctive use, but within the area of influence of the 
groundwater pump. In developing the model, the volume 
of groundwater is assumed constant, so that the depth of 
groundwater pumped (P) from the area (1 — U) A, outside 
the conjunctive use area is equal to groundwater recharge 
in that area as described by Egq.(5). The salinity of 
groundwater outside the conjunctive use area, C,, ,can be 
described by Eq. (6), which is an expression for ground- 
water salinity where drainage water is pumped but there 
is nO conjunctive use on the soil surface above. Substitut- 
ing Eq. (6) for C,,, in the abbreviated form ae~>'+¢, and 
Eq. (5) for P into Eq. (9) yields: 


dC,, _ +R) LF U(W, +S)C, +RC,) 


dt 


H-60-U 


Integrating Eq. (10) and using the initial condition C,,= Cy at t = 0 yields: 


[(W, + R) LF +S] + ((W,+S)C,+RC,)t (11) 


+ 


The rate of increase in groundwater salinity in the 
conjunctive use system in Fig. 4, dC,,/dt, can be ex- 
pressed by adding the salts entering the area from adja- 
cent pumped groundwater, (1—U)A,, to that from rain- 
fall plus irrigation supply water applied and channel seep- 
age in the area of conjunctive use, and dividing by the 
volume of groundwater under consideration. That is: 


dC,, _ PC,,(1—U) A, + (W, +S)C, + RC) UA, 
HOUA, 


(9) 


Equation (11) is an expression for groundwater salinity as 
a function of time under spatially constricted conjunctive 
use. 

For full irrigation the depth of irrigation supply water, 
W,, applied in the conjunctive use area depends on the 
volume of groundwater pumped from outside the con- 
junctive use area. Thus W, in Eq. (11) is a function of U. 
Such an expression for W, can be derived from the mass 
conservation of water and estimated crop water use for 
this system. For full irrigation, maximum crop evapora- 
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tion, ET,, = 0.85 E,, can be expressed as a function of 
channel supply water (W,), groundwater pumped from 
the area of conjunctive use (J - LF + S), the depth of water 
lost to leaching (J - LF), in the area of conjunctive use, 
and the groundwater pumped from outside the area of 
conjunctive use, as expressed by Eqs. (5) and (8). The 
depth of water use can then be calculated by summing 
water volumes and dividing by the area of conjunctive use 
to yield: 


UA(W,+R+4+I1- 


0.85E,- LF 
1—LF 
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been spatial constriction of conjunctive water use will 
also reverse the trend of increasing groundwater salinity 
caused by spatial constriction of conjunctive use. 
Problems introduced by spatial constriction of con- 
junctive water use are likely to be greatest where ground- 
water salinity is marginal, or where surface water supplies 
are inadequate, and there is therefore greater incentive to 


ET, = 0.85E, = 


Rearranging terms to express W, in terms of U yields: 


W, = 085E, — 


+s|-R-s. 


(13) 


Equation (13) can be used to calculate the applied depth 
of channel supply water required to be mixed with 
groundwater for maximum crop evapotranspiration in 
the area of conjunctive use in Eq. (11). Equation (8) can be 
used to estimate W, for Eq. (11). Partial irrigation is rep- 
resented by use of a crop factor of 0.75 instead of 1 for 
Fig. 2. The predictions of Eq. (11) are then plotted in 
Fig. 5 with parameter values typical of the Shepparton 
Irrigation Region in the Murray Basin, Australia. Fig- 
ure 5 illustrates that very rapid degradation is possible 
with low values of U, demonstrating the importance of 
distributing pumped groundwater over the entire area of 
influence of the groundwater pump. 

The rate of groundwater degradation from spatially 
constricted conjunctive water use is higher with smaller 
depth of groundwater in the system, and with greater 
values of the leaching fraction. The rate of degradation is 
also greater with higher levels of initial groundwater 
salinity (see Eq. (11) and Fig. 5). Partial irrigation reduces 
the impact of spatially constricted conjunctive water use 
(Fig. 5) because of the reduction in leaching fluxes result- 
ing from the reduction in J. Comparison of partial and full 
irrigation in Eq. (11) and Fig. 5 for the same value of LF 
underestimates the effect of partial irrigation since a re- 
duction in J could also be expected to reduce LF. 

Pumped groundwater salinity may increase excessive- 
ly after some time due to the effect of spatial constriction 
of the area of conjunctive use. The increase could be coun- 
teracted by installation of another groundwater pump in 
the area of groundwater outside the area of conjunctive 
use, where groundwater salinity will be lower. Operation 
of the pump subsequently installed in the area of low 
salinity groundwater will then result in groundwater salt 
transport out of the area of higher salinity groundwater 
to the area around the pump. Problems with this strategy 
include the cost of exploration for additional suitable 
pump sites, the potential lack of available pumpable 
aquifers for pump relocation, and the lack of on-farm 
management control if the additional pump site exists 
beyond the farm boundary. Application of pumped 
groundwater over the whole of the area of influence of the 
groundwater pump after some period where there has 


(12) 


pump groundwater to augment surface supplies. When 
groundwater is used over the whole of the area of influ- 


ence of the groundwater pump, the depth of groundwater 
for conjunctive use equals the depth of recharge, which is 


usually a small component of the surface water budget. 
Conversely, if a substantial proportion of crop water use 
is supplied by pumped groundwater, a spatial constric- 
tion of conjunctive water is implied, and excessive degra- 
dation will occur. 


Groundwater mixing 
The salinity of pumped groundwater varies over time 


because of spatial variability in groundwater salinity and 
variations in travel time from the point of recharge to the 
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Fig. 5. Groundwater degradation from spatially constricted con- 
junctive water use. Plot of Eq. (11) with W, and W, described by 
Eqs. (8) and (13) respectively. Parameters nominated are: C, = 
0.1 dS/m, R=0.45 m/y, C,= 0.008 dS/m, = 0.04 m/y, 0.4. 
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pump screen. This variability will reduce over time as 
groundwaters mix. When a pump is installed in ground- 
water surrounded by better quality groundwater, pumped 
groundwater salinity will improve over time through mix- 
ing, and vice versa. When groundwater salinity is being 
assessed for conjunctive use, the total storage within the 
area of influence of the groundwater pump should be 
considered, so that the effects of mixing on groundwater 
degradation can be assessed. 

Spatial variations in groundwater salinity will not affect 
the underlying processes of groundwater degradation 
already discussed, but it is possible that the resulting tem- 
poral variation in pumped groundwater salinity will 
mask these processes when field data is assessed. In some 
cases spatial variability in salinity may be so great that 
specific management strategies may be necessary to 
achieve optimal salinity control. For example, reduction 
in the pumped water depth through intermittent pump- 
ing, or reduction in recharge through partial irrigation, 
will increase the time taken for upper groundwaters to 
reach the pump screen. Such action could be beneficial in 
areas with higher groundwater salinity in the upper stra- 
tigraphy. However, intermittent pumping may result in 
some loss of leaching and perhaps productivity, as out- 
lined by Prendergast et al. (1994). Where underlying 
groundwater is of poor quality a reduction in the rate of 
groundwater deterioration could be achieved by reducing 
the pumping rate, which reduces well drawdown and 
therefore the hydraulic gradient from deeper ground- 
water to the conjunctive use system. 


Conjunctive water use incorporating limited export 


The effect on groundwater salinity of off-site disposal of 
all drainage water was described in Fig. 3. Long term 
sustainability of conjunctive use is also affected if salt 
export is reduced and made equal to the mass of salt that 
enters the irrigation area with irrigation water and rain- 
fall; this will be referred to as export for salt balance. 
Prendergast et al. (1994) suggest that export for salt bal- 
ance has little effect on crop yield in the conjunctive use 
system in the short term, because this export usually con- 
stitutes only a small fraction of the groundwater pumped. 

Equation (2) describes the effect of rainfall and irriga- 
tion salt inputs on groundwater salinity under conjunc- 
tive use; if export is equal to surface salt inputs for this 
case groundwater salinity remains constant. The rate of 
degradation from salt inputs from irrigation and rainfall 
depends on the depth (H) of the conjunctive use system 
(Fig. 2). Thus export for salt balance results in greater 
reduction in groundwater degradation when shallow 
aquifers are utilised (small H), and has little impact when 
A is large. 

The effect of export for salt balance on the groundwa- 
ter system described by Eq. (4) and Fig. 3 with positive G 
(which results in a rate of degradation less than that 
caused by irrigation water salt inputs) is to reduce 
groundwater salinity over time. For a negative value of G, 
referred to in the section “Upward leakage into the con- 
junctive use system”, export for salt balance will reduce 
the rate of groundwater deterioration. 


--- With Export 
—— Without Export 


Groundwater Salinity , C,, (dS/m) 


4 
T 


100 150 200 


Time, t (Years) 


Fig. 6. The impact of salt export (made equal to the rate of surface 
salt additions) on spatially constricted conjunctive water use. 
Broken lines represent the conjunctive use system with export 


The effect of off-site export for salt balance on the 
groundwater system described by Eq. (11) (spatially con- 
stricted conjunctive water use) can be assessed by sub- 
tracting the impact of export; the impact of export on C,, 
can be quantified as the annual salt inputs from irrigation 
water and rainfall (ie. (W,C,+RC,)(1—U)+(W,C, + 


RC,) U) divided by the water depth H 6 and multiplied by 
t. Equation (11) has been plotted in Fig. 6 as a broken line, 
after subtraction of this expression for export from the 
right hand side of Eq. (11). 

Export for salt balance has the effect of reducing 
groundwater degradation caused by spatially constricted 
conjunctive water use to zero in the long term, after initial 
rapid degradation. The time to reach this equilibrium 
depends on the depth (H) of the conjunctive use system. 
Equilibrium in the conjunctive use area occurs when the 
groundwater salinity within the area of influence of the 
groundwater pump, but outside the area of conjunctive 
use, comes into equilibrium with the salinity of ground- 
water recharge there (i.e. approaches steady state). Initial 
degradation from spatial constriction of conjunctive use 
depends on the difference between the depth weighted 
groundwater recharge salinity and Cy, and is therefore 
greater at higher levels of Cy, irrespective of export for 
salt balance. The difference between the final steady state 
groundwater salinity and the initial groundwater salinity 
in Fig. 6 is also greater for higher initial groundwater 
salinity and smaller values of U. 

The rate of groundwater degradation depends on the 
depth of groundwater in the conjunctive use system 
(Eqs. (2), (4), (6) and (11)), spatial distribution of irrigation 
water and groundwater (Eq. (11)), the spatial relationship 
between the regional groundwater and the conjunctive 
use system and regional groundwater flows which deter- 
mine G (Eq. (4)), groundwater mixing, water management 
(including regional redistribution of groundwater), dis- 
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posal, and the depth of non-leaching recharge (e.g. S$ in 
Eqs. (4), (6) and (11)). If groundwater degradation is to be 
minimised, management practices which influence these 
factors must be optimised; disposal for salt balance is 
insufficient to offset inappropriate groundwater manage- 
ment. 

Problems likely to be encountered with export for salt 
balance include the downstream environmental costs 
when a river receives the effluent, the cost of drainage 
infrastructure, and difficulties in accurately calculating 
the salt mass for export from a groundwater pump, be- 
cause this mass depends on the pump’s area of influence, 
which may be difficult to determine. Difficulties in deter- 
mining a pump’s area of influence arise because of hetero- 
geneity on local stratigraphy, and temporal variations in 
pumping and groundwater recharge rates and because of 
elements of confined aquifer behaviour; in locally con- 
fined aquifers, under constant pumping the area of influ- 
ence expands over time until the recharge zone is reached. 


Model testing and groundwater degradation 
in the Shepparton Irrigation Region 


Since the advent of high watertables in the Shepparton 
region late in the 1970 decade, there has been a steady 
expansion in the area affected by salinity. A 640 ha 
project area was established at Tongala in 1981 to moni- 
tor the effects of conjunctive use of groundwater on pro- 
ductivity and salinity (Heuperman 1988). In the project 
area 90% of supplied irrigation water is applied to irrigat- 
ed perennial pastures that support dairy production, 
which dominates irrigation water use in the Shepparton 
region. Conjunctive use of groundwater has been shown 
to result in considerable reduction in soil salinity in the 
area (Heuperman 1986). 

Average groundwater degradation in the Tongala 
project area from irrigation water salt inputs and rainfall 
would be expected to be approximately 0.012 dS m~! per 
year from Eq. (2), if a lower confining layer at 10 m is 
considered to be base of the groundwater streamlines. 
Pumped groundwater salinity has increased by approxi- 
mately 0.07 dS m~' per year over eight years of measure- 
ment (Prendergast and Heuperman 1988). Second order 
finite element modelling indicated that approximately 
40% of the pumped groundwater volume flowed into 
the project area from outside the project boundary 
(H. Kleindienst pers. comm.). Difficulties encountered in 
the modelling included the allocation of appropriate 
boundary conditions and the obtaining of sufficient de- 
scriptive data. Equation (11) can be used to estimate 
groundwater degradation which could be expected if the 
lower confining layer is considered to be the base of the 
streamlines. Suitable parameters are C, = 0.1 dS/m, 
C,=0.008 dS/m, R=045m, 6=04, H=10m, 
C, = 3 dS/m, LF = 0.1, W, = 0.6 m, W, = 0.5 m, S = 0.04, 
and U = 0.6. Evaluation of Eq. (11) using these values 
yields an average groundwater degradation of 0.1 dS/m/ 
yr. It is difficult to assess the accuracy of this prediction 
given the long time scale of the processes involved, as 
illustrated in the previous section, and the short time scale 
of the measured data. 


Discussion 
Groundwater pump management 


Aspects of groundwater pump installation and manage- 
ment can also have a substantial impact on groundwater 
degradation. The extent of the area of influence of the 
pump is not always easily determined because of hetero- 
geneity, and many extend across farm boundaries. Farm 
boundaries can restrict the area of conjunctive use, but 
not A,, so that reduced U and increased groundwater 
degradation can result (Eq. (11)). The depth of groundwa- 
ter flow per unit area entering the conjunctive use area 
from outside is likely to be less when a network of pumps 
operates over a large area. A large area of conjunctive use 
means that the ratio between the perimeter of the area of 
influence and the area of influence is less than for an 
isolated pump. 

For isolated pumps, reduction in groundwater flows 
into the conjunctive use system, and therefore a reduction 
in groundwater degradation (and an increase in U) could 
be achieved by intermittent pumping. This reduces the 
time that lower hydraulic potentials exist in the conjunc- 
tive use system and therefore reduces the depth of 
groundwater entering the system from outside the con- 
junctive use area. Intermittent pumping will, however, 
have the effect of restricting leaching which may result in 
increases in rootzone salinity, despite smaller surface ap- 
plied salt loads (Prendergast et al. 1994). 


Catchment management 


Optimal salinity management of conjunctively used water 
can only be achieved through effective groundwater 
catchment management. The conjunctive use system is 
invariably part of a regional groundwater flow system, so 
that the boundary conditions for the conjunctive use sys- 
tem are affected by regional groundwater. Lateral flow 
can add or deplete salts from the pumped aquifer and 
therefore change groundwater salinity. Where regional 
scale recharge/discharge processes occur (Hubbert 1940; 
Fetter 1988), surface redistribution of pumped ground- 
water can enable some groundwater leakage to be main- 
tained in the conjunctive use system over the whole irriga- 
tion region, if net groundwater outflows occur on 
average. In this way redistribution reduces excessive 
groundwater degradation where discharge from the re- 
gional groundwater to shallow aquifers occurs. Ground- 
water salinity in the conjunctive use system is then de- 
scribed by Eq. (4). Also, conjunctive water use relies on 
crop threshold salt tolerance to facilitate surface disper- 
sion of groundwater without resulting in crop yield loss; 
groundwater redistribution can enable utilisation of crop 
threshold salinity tolerance over as large an area as possi- 
ble (Prendergast et al. 1994). 

Conjunctive use without groundwater export at the 
surface will probably result in higher regional soil and 
groundwater salinities which will approach a steady state 
only after thousands of years. The higher regional equilib- 
rium will result in greater natural discharge of salts from 
the landscape in comparison to salt discharge at earlier 
times (as illustrated at a smaller scale by Eq. (4)). Greater 
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salt discharge will occur through both regional ground- 
water outflows and surface runoff to rivers. It is impossi- 
ble to assess what the increases in regional salt storage are 
likely to be, and whether changes in agricultural practices 
and salt discharge from the landscape will be sufficient to 
offset the effects of higher salinity over these long time 
scales. 

The implications of conjunctive water use for the long 
term could also be evaluated by assuming there is no 
increase in salt discharge from an irrigation region over 
time. In the Shepparton region, for example, there is net 
natural flow out of the groundwater system, although 
local aquifers can experience either recharge to or dis- 
charge from deeper aquifers. There is potential for in- 
creases in salt storage within both the shallow and deep 
aquifer systems. The effect of utilising this storage to re- 
duce salt accumulation on the soil surface, on a regional 
scale, can be illustrated by a simple example which ap- 
proximates conditions in the Shepparton region. If the 
irrigation area comprises 0.5 M ha and has average val- 
ues of H = 100 m, 6 = 0.4 and C, = 3 dS/m, the mass of 
salt storage equals (0.5 x 10° x 10*) m? x 100 m x 0.4 x 
(3 x 0.6/1000) = 3.6x 10° t, by using the generally ac- 
cepted approximation 1 dS/m = 600 mg/l. If the irriga- 
tion area has a net salt addition of 50,000 t/y (as is ap- 
proximately the case with the Shepparton region), the 
resulting percentage increase in regional salt storage in 
100 years is approximately (5 x 10* x 100/3.6 x 108) x 100, 
or 1.4%. Conversely, export of 50,000 t/y of salt from the 
region will reduce regional salt storage by 1.4% after 
100 years, in comparison to having no export available. 
The effectiveness of conjunctive water use in this case will, 
therefore, ultimately be determined by integrated land 
and water management, which governs the distribution of 
salts in the landscape and the groundwater systems, and 
can prevent excessive salt accumulation in the rootzone. 


Conclusions 


Under optimal land and water management, salinity at 
the base of the rootzone in a conjunctive use system will 
ultimately be determined by net groundwater flows out of 
an irrigation area, if there is no export of salt (Eq. (4) and 
Fig. 3). Because irrigation regions usually have some net 
natural groundwater outflows, groundwater degradation 
to higher equilibrium salinity occurs slowly over many 
centuries, if good management of the groundwater catch- 
ment is maintained. Good salinity management of the 
groundwater catchment includes strategies such as: 

(a) Distribution of groundwater over the entire area of 
influence of a groundwater pump. Pumping groundwater 
to obtain a water depth greater than recharge, to augment 
irrigation, will induce faster degradation and can be con- 
trary to good salinity management (Fig. 5). 

(b) Utilisation of maximum depth of stratigraphy for 
salt storage through deep installation of groundwater 
pumps (Fig. 1). Increasing the depth of groundwater for 
salt storage reduces the detrimental effects of salt inputs 
from rainfall and irrigation (Fig. 1), discharge from the 
regional groundwater, groundwater mixing and spatial 


constriction of conjunctive use (Fig. 5). Deep installation 
must be weighed up against adverse affects of poorer con- 
trol of waterlogging and of groundwater mixing if the 
deeper stratigraphy contains layers of poor quality 
groundwater. 

(c) Reduction of salt inputs from irrigation water, 
achieved through partial irrigation and reduction in irri- 
gation intensity (Fig. 1). Partial irrigation increases the 
area which can be irrigated. Partial irrigation also re- 
duces the impact of spatial constriction of conjunctive 
water use (Fig. 5) by reducing leaching fluxes both 
through a reduction in total applied water and a probable 
reduction in the leaching fraction. 

(d) Surface redistribution of pumped groundwater 
to minimise the effects of regional scale groundwater 
recharge/discharge conditions caused by heterogeneity in 
topography and hydrogeology. Installed pumps in a con- 
junctive use area can provide leaching for irrigation in 
distant areas, which results in high negative G in the con- 
junctive use area, and therefore excessive degradation, if 
there is no redistribution. 

(e) Development of specific groundwater pump man- 
agement strategies, such as intermittent pumping, in areas 
where groundwater mixing is likely to cause excessive 
degradation, or where only isolated pumps can be in- 
stalled. These strategies must be developed in the context 
of local agricultural practices, climate and stratigraphy. 

Artificial means of exporting saline groundwater out 
of the irrigated area also reduces groundwater degrada- 
tion. However, export from each installed pump, equal to 
salt additions at the surface, can be costly in terms of 
infrastructure and downstream environmental impact, 
and may only significantly affect groundwater degrada- 
tion where just shallow aquifers can be used to control 
salinity. Where groundwater salinity is too high for con- 
junctive use, groundwater export will be necessary if max- 
imum productivity is to be maintained. 

Data to adequately describe groundwaters in detail 
for a specific site are difficult and costly to obtain. This is 
especially the case where high spatial variability in strati- 
graphy exists, such as is the case in the fluviolacustrine 
deposits of the Murray Basin. Because of the long time 
scale of the processes involved, at least several decades of 
intensive in-situ measurements will be required for better 
model validation. Nevertheless, the basic principles of 
mass conservation of salt and water can be applied as 
illustrated to provide management guidelines to enable 
sustainable conjunctive water use. 
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Abstract. A kinematic-wave model is developed for simu- 
lating the movement of soil moisture in unsaturated soils 
with plants. The model involves three free boundaries. 
Analytical solutions are derived when the plant roots are 
assumed to extract moisture at a constant rate and the 
upstream boundary condition is independent of time. 
Numerical solutions are the only resort when the mois- 
ture extraction and the upstream boundary condition 
both depend on time. 


Much of the mathematical treatment of flow in unsatu- 
rated porous media has dealt with capillary-induced flow 
(Smith 1983). However, there exists a multitude of cases 
where gravity dominates vertical movement of soil mois- 
ture. Some examples include: drainage following infiltra- 
tion, the water percolation deeper into the soil, the verti- 
cal movement of moisture in relatively porous soils when 
rainfall or surface fluxes are typically of the order of, or 
less than, the soil-saturated hydraulic-conductivity, to 
name but a few. For treatment in such cases, the 
kinematic-wave theory is simple yet reasonably accurate. 

Although Sisson et al. (1980) applied the kinematic- 
wave theory to internal drainage, Smith (1983) was prob- 
ably the first to apply the theory to develop a complete 
kinematic-wave model for soil moisture movement. 
Charbeneau (1984) extended Smith’s work to solute 
transport, and Charbeneau et al. (1989) to multiple solute 
transport. In a series of papers, Germann and coworkers 
(Germann 1985; Germann and Beven 1985, 1986; Ger- 
man et al. 1987) extended the application of the theory to 
infiltration and drainage into and from soil macropores, 
as well as to microbial transport. Yamada and Kobayashi 
(1988) discussed the kinematic wave characteristics of 
vertical infiltration and soil moisture, with the aid of field 
observations on tracers. They concluded that the vertical 
infiltration of soil moisture had the characteristics of 
kinematic waves. 


Correspondence to: V. P. Singh 


In this study, a kinematic-wave model is formulated 
for soil-moisture movement with plant-root extraction. 
Analytical solutions are derived under the condition that 
the moisture extraction by plants and the surface flux are 
both constant in time. This condition is severe, but leads 
to intriguing and useful results. For more realistic con- 
ditions numerical solutions are the only resort, these 
together with field verification, will be reported in the 
near future. 


Kinematic-wave model 


Formulation 


For vertical unsaturated flow with plant-root extraction 
and with the neglect of the flow of air, the governing 
equations are the law of conservation of mass and a flux 
law. The one-dimensional conservation of mass equation 
or the continuity equation can be expressed as 


(1) 


where @ is the volumetric moisture content (dimension- 
less, volume of water/volume of soil), q is the vertical flux 
of soil moisture (L/T), e(z, t) is the rate of soil moisture 
loss due to plant-root extraction at time t (1/T), z is the 
depth (L) measured downward from the ground surface, 
t is the time (T) that water has stayed at location 
t=t—w/(z), w(z) is the time history of the moisture ad- 
vance front, and t is time (T). At a fixed t, e(t) is essen- 
tially the rate of plant-root extraction of soil moisture on 
a unit depth basis defined as 


= e(z, T) (2) 
n 


where n is the soil porosity, and r(t) is the soil-moisture 
extraction rate per unit depth at time ¢ (1/T). It may, 
however, be reasonable to assume e(z, t) = e(t). This as- 
sumption implies spatially uniform extraction, that may 
be true if the plant roots are uniformly developed in the 
soil. 
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The flux law for a kinematic-wave model can be ex- 
pressed as 


= (3) 


where K(@) is the unsaturated hydraulic conductivity 
(L/T) as a function of the soil-moisture content. A vari- 
ety of expressions have been proposed for K (0) (Hsu and 
Liu 1990). Of all the expressions, the best-known perhaps 
is the Brooks-Corey relation (Brooks and Corey 1964) 
which can be written as 


K(0)=K, ( (4) 
where K, is the saturated hydraulic conductivity (L/T), 
0, is the residual value of @ below which water cannot be 
extracted (or moved) by capillary forces, @, is the saturat- 
ed water content = porosity = n, and a is a parameter 
typically between 3 and 4, and is related to the pore-size 
distribution index. 

In Eg. (1), t is unknown and a relation is needed 
whereby t can be determined. With such a relation, the 
formulation of the kinematic-wave model will be com- 
plete. To that end, when water is applied to the soil at its 
surface, it infiltrates into the soil and advances vertically 
downward, increasing the moisture content. If the soil is 
initially dry, the advance front, also called the shock 
front, will define the interface between moist and dry 
parts of the soil matrix. Let z = s(t) or t = w(z) be the 
time history of that advancing front; this time history 
gives the space-time coordinates of the shock front. The 
front is a free boundary and has to be determined along 
with the solution. An equation for z = s(t) or t = w(z), 
therefore, has to be formulated. This can be accom- 
plished by observing that 


0) K(@ 
u(6) = u(z,) =) = (5) 


where u(@) is the average velocity with which @ moves. 
Thus, the free-boundary equation is obtained by replac- 
ing z by s(t) or t by w(z) in u(z, t) as 


ds(t) _ K(O(s(), 0) 

or 
dw(z) _ -1_ 9(z, w(z)) 

dz ix [u(z, w(z))] K (0(z, w(z))) 
Equation (6) is valid for 0(z, t) > 0. Equations (6) and (7) 
require the advance front to move with the speed immedi- 
ately behind the front. 

The kinematic-wave model formulation consists of a 
partial differential Eq. (1), an algebraic Eq. (3), and an 
ordinary differential Eq. (6) or (7) with two unknowns 
0(z, t) and s(t) or w(t). In order to derive solutions of these 
equations, the following can be assumed: 


A(z,0)=f(z), z>0, (8) 
0(0,t)=g(t), O<t<T (9) 
=0,, t2>T 


s(0)=0 (6) 


w(0)=0. (7) 


where T is the duration for which moisture is applied at 
the upstream boundary. Note that g(0) = f(0) = 0). 
Equations (1) and (3) can be combined to produce 


00 dK(0) 00 


For notational simplicity, let 


dK(0) dz 
11 
dé 00/ Oz ot constant 
which is referred to as the mobility of water in soil by 
Irmay (1956); this is analogous to celerity, commonly 
used in open channel flows. 


Solution 


Equation (10) is satisfied in the domain bounded by the 
positive t axis and the curve t = w(z) or z = s(t). Solution 
of Eqs. (10) and (6) can be derived by using the method of 
characteristics. Accordingly, one can choose o as the 
parameter on the t-axis and t as the parameter along 
the characteristics. The characteristics originate from the 
t-axis on the segment 0 < t < T. The solution domain can 
be divided into two subdomains D, and D,, with the 
characteristics z(t, T) serving as the dividing line (called 
the bounding characteristic), as shown in Fig. 1. The 
characteristic curves z(t, 0(t,o) passing through the 
points (0, o, g(a)) in the (z, t, 0) space, satisfy 


d0(t, o) 


dz(t, 
dt 


= M(t,a), z(o,c0)=0. (13) 


P=(z,,t,) 


Z =2 (t,o) 
(§ (6), 


0 


Fig. 1. Characteristic solution domain for unsaturated flow with 
plant-root extraction 
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z(t, o) is the position at time t of the soil-moisture content 
@ which was at z = 0 and at t = o; O(t, ) is the soil-mois- 
ture content at time t. The plant-root extraction rate de- 
pends upon the nature of vegetation and availability of 
soil moisture. If e(t) is constant, which is not entirely 
unreasonable, especially when the soii-moisture is not 
limiting plant root extraction, then e(t) = e. 

Depending upon the nature of the function g(t), three 
cases can be distinguished: (1) g(t) = constant, (2) dg(t)/dt 
<0, and (3) dg(t)/dt >0. To describe the variation of soil 
moisture in time at a fixed z or in space at a fixed t, one 
may determine 00/0z and 06/0t. 


Solution for domain D, 


This domain is bounded by 0<t<T, z(t, T), and 
z = s(t). Solution of Eq.(12) is 


A(t, 0) = g(o) —(t—o)e. (14) 


Substitution of Eq. (14) in Eq. (13), and then integration 
yields 


z(t,o) = M@(o)—(t—9) e)dt. (15) 
With the use of Eq. (4), 


K,a (16) 
(0, — 95) \ 8, — 9% 


Substitution of Eq. (14) in Eq. (16) yields 


aK, 
(8, — | 0, — 9% 


With the use of Eq. (17), Eq. (15) leads to 

(18) 


Equations (14) and (18) constitute the characteristic solu- 
tion for domain D, . By eliminating o between these equa- 
tions, 8 can be expressed as a function of z and t. To do 
that, z(t, o) must be, for fixed t, either an increasing func- 
tion of o or a decreasing function of o. Differentiating 
Eq. (18) with respect to o yields 


—§ 
a (a-1)/a 
g(a) 
0, K, 


e(0(t, 6) — (20) 


It is seen from Eq. (20) or (19) that g(a) = dg/do < Oisa 
sufficient condition for z,(t,o)<0. The condition 


M(6)= 


o) = 


(17) 


or 
aK, 
(0, 


z,(t, 0) = 
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g(a) <0 includes the case of principal interest to us, 
namely, g(a) = constant. 


Case 1. g(t) = constant: If g(o) = 0,, then Eq. (18) be- 


comes 

K, 0,—9% |° 

By coupling Eqs. (41) and (21), @ can be expressed as a 
function of z and t as 


a l/a 
0(z) = 05+ 4) 


It is seen that in domain D, , 0 is a function of z and does 
not depend on t. The soil moisture profile is nonuniform 
but steady-state. In Eq. (22a), there are two terms within 
braces. The second term is due to water uptake and the 
first term is due to the boundary flux. If e=0, then 
0(z)=6,. The nonuniformity of soil moisture is a direct 
consequence of the plant-root extraction. If e were zero, 
then 6 would be independent of z. Thus, the alteration of 
the soil moisture profile is caused by plant roots. Of 
course, the degree of alteration depends on the value of e. 
On a short-time scale, the value of e is relatively small in 
field. Equation (22a) can be written with use of Eq. (4) as 


a l/a 
6(z) = 0, + (0, — 


In many practical cases, e z/K (8,) is much less than 1, and 
attenuation of the 6-profile, as a result, will be small. 


When Eq. (22a) is inserted into Eq. (4), the soil-moisture 
flux, with the aid of Eq. (3), is obtained: 


0, — 0,1" 
K (0) = q(z) = K, (23) 


z(t, o)= 


(22a) 


Similarly, the average velocity of the soil-moisture move- 
ment is obtained by substituting Eq. (22) in Eq. (5), 


|? ez 


Similarly, substitution of Eq. (22) in Eq. (16) yields the 
mobility of water: 


s 


u(z) = 


(24) 


The behavior of the soil moisture in this case can be fully 


described. From Eq. (22) at a fixed z, z > 0, 


00 
26 
ot 


and at any t,0<t<T, 


(1-a)/a 
— 
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Equation (27) shows that at any depth within D,, 6 re- 
mains constant in time and at any time it decreases with 
increasing z from the ground surface. The gradient of 0 
increases with increasing z. This corresponds to the case 
of a wetting event. 


Case 2. dg(t)/dt < 0: In this case, 00/0z and 00/Ot are ob- 
tained from Eqs. (14) and (18) by noting that 


dt 


To that end, with g’= g(o) = 0g(a)/do, 


(28) 


e (0, — 05) 


(a-1)/a 


g — 9% ez (1 -aj/a 
(31) 
Therefore, by dividing Eq. (29) by Eq. (30), 
32 
00 + e)(9, — 65) (62) 


{t should be noted that in this characteristic method of 
solution, ¢ was chosen as the parameter along character- 
istics. Therefore, in Eq. (18), z is the dependent variable, 
and o is the independent variable. In order to obtain 
00/Ot, one can choose t as the dependent variable and z as 
the parameter along characteristics. The procedure, how- 
ever, remains the same in both cases. Employment of t as 
the dependent variable results in 


and 


ot g a— a ez 

(32b) 
By dividing Eq. (32a) by Eq. (32b), 


The variation of 6 with z for a fixed t can be analyzed 
from Eq. (32). Ife < g’(a) then 00/dz > 0; that means that 
6 increases with increasing z. On the other hand, if 
e=g'(c), then @ is independent of z. This is also the case 
when z > If e > g'(a) then de- 
creases with z but can also have a complex variation. At 
a fixed z, the variation of o can be analyzed with the aid 


of Eq. (33). If g(a) =0, then 6 is independent of t. If 
g(a) <0, then 6 decreases with time. If e = g'(c), then 6 is 
a decreasing function. At z = K,((g—0,)/(0,—0,))*/e, 
then 


(34) 


00 e(g'+ e) 


Equation (34) shows that 0 becomes independent of t 
whereas Eq. (35) shows that as a function of z, @ decreases 
with z. 


Case 3. dg(t)/dt > 0: This can be analyzed from Eqs. (32) 
to (35). However, Eq. (20) reveals that when g’(c) > 0, 
z,(t, ¢) > 0 and this case may entail shock formation. In 
that case, the analytical treatment becomes unwieldly and 
is not discussed in this paper. 


Determination of advance front (free boundary 


of domain D, ) 


In order to determine the shock front t = w(z), it is conve- 
nient to express it in terms of the parameter o introduced 
above. If one considers the characteristic curve z(t, o) in 
the (z, t) plane, then it will intersect the free boundary at 
the point denoted by (&(c), 4(c)) as shown in Fig. 1. 
Therefore, z = €(c), t = n(c) is the parametric representa- 
tion of the free boundary in terms of the parameter o. For 
the case g(t) = 0, =constant, and e(t) = e = constant, 
determination of the free boundary is relatively simple. 
From Eq. (22), with t replaced by (a), 


2(n(o), = (36) 
K, 8, — 8, — (n(o) — o)e — 8, 
From Eggs. (6), (14) and (4), and noting that 6 does not fall 
below 


ds(t) _ K(0(E(o), n(o)) 
dt (a), n(o)) 


K, 6, —(n(o) — o)e — 
0, — 9, | 


(37) 


Equations (36) and (37) can be solved for &(c) and n(a). To 
that end, Eq. (36) can be differentiated to yield 


ak, — — a)e— 
(0, — 0, — 


O(a) = 


(38) 


On eliminating ¢’(a) between Eqs. (24) and (25), 


(0) =0 (39) 


which gives 


(40) 


a-— 


= 
00 
00 | 
K, 
Ang 
| 
| 
| 
q 


Similarly, eliminating n'(c) between Eqs. (37) and (38) 
leads to 


which gives 
Ks 
9 eo j 


Equations (40) and (42) constitute the parametric repre- 
sentation of the free boundary. By eliminating o between 
Egs. (46) and (42), and noting that € corresponds to z and 
n to t. In terms of q, 


_ K, do 1/a et a 
where qo= q(@,). Equation (43) or (44) expresses the ad- 


vance front. It is seen from Eq. (43) that the coordinates 


of the point to which the free boundary can extend is 
given by 


_K, 


To complete the solution for domain D, , the intersec- 
tion of the bounding characteristic z(t, T) and the free 
boundary z = s(t) is to be determined. Let that point be 
P =(z,, ts). This point can be determined explicitly by 
equating Eq. (21) with o = T, t = ty and z = z, to Eq. (43) 
with z = z, and t = t,; that results in 


(46) 


(47) 
For the case when g(t) is not constant, we use the same 
parametric representation as before. Therefore, analogous 

to Eq. (36), 

E(o)= {Ig (e)—9}*— [g(e) —(n(e) — 0) e— 
(48) 


and analogous to Eq. (37), 
n'(o) (0, — 
Differentiating Eq. (48) and then substituting the differen- 
tial in Eq. (49) yield 

a g (2) [g(o) — 

— (n(o) — 9) e — 


nO=0. 


(49) 
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Equation (50) cannot, in general, be solved in terms of 
simple forms of g(t), an analytical solution for y(c) is 
tractable. With substitution of this solution in Eq. (48), 
€(a) is obtained. If, however, g' (a) = 0, Eq. (50) reduces to 
Eq. (39). 


Another method to determine the advance front 


The velocity of the front can be expressed as 


uls(,)= 


(51) 


Representing 0, by 0, and 0, by 0, and taking advantage 
of Eqs. (22) and (23), Eq. (46) can be written for the case 


g(t) = constant, as 


as()__ Ki" 
dt 0,— 4 


which, when integrated, leads to 


K, q \'* te 


which is the same as Eq. (44). 


(52) 


z=s(t)=— 


Solution for domain D, 


Domain D, is bounded by T <t, z = z(t, T), z = s(t) be- 
yond the point p = (zx, ts). In this case, 0(0, t) = 0, for 
t>T; Let g(t, be a continuous function of t 
coinciding with g(t)on0<t< 
as shown in Fig. 2. Then, it is seen from Eqs. (14) and (18) 
that the characteristics z = z(t, 0, 6), T<o < T +6, cover 
the region above z = z(t, T ) completely. 


O(t, o, — 9, |? 
0,— 9 | 


z(t, o, = 


~» Fig 2. Function g(t, 


T+€E 


| 

| 

| | 

| 

| 

| 

ty = | | 

{ a-—1 

| 
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By letting «+0 and o > T, one obtains 


| 
This equation expresses the soil-moisture 6 as a function 
of z and t in domain D,; 6 varies from 0, to 6). The 
reduction in @ is caused by the plant-root extraction, and 
hence the water uptake is the dominant factor for varia- 
tion of water content. This can be easily shown for a 


special case when a is an integer. To that end, Eq. (54b) 
can be written as 


2(0, t) = (; + - i}. (55a) 


With use of the bionomial theorem, this can be written as 
_ K, (0-95 \" (a\ (et -—T)\! 


It is seen that each term in the series is affected by the 
plant-root extraction e raised to the power (a—i—1). 


Assume that only the first term of the series is most signif- 
icant. Then, 


ae (55c) 


Thus, the reduction in @ is primarily caused by water 
uptake. 


The characteristics in domain D, are given by 


_ K, 0;— 0;-—e(t— T)-—0,\" 

e ( 0,— (6) 
where 6; is the initial soil-moisture 0)<0;<0,, that is 
carried out by the i-th characteristic. 


Determination of free boundaries for domain D 


Domain D, has two free boundaries FB2 and FB3 as 
shown in Fig. 3. Along FB3 that starts at t = T and x = 0, 
0(z, t) = 0); thus, this is the locus defined by 0, . This is the 
time history of the drying front as it advances from x = 0. 
From Eq. (14), we have 


T<o<Tte. (57) 
Thus, from Eq. (18), we obtain 
z(t, o, 8) = : (58) 
Letting 0, o > T, we get the free boundary FB3 


Between z=0 and Eq. (59), 0(z, t) = 0. 


FB3 
Z P= (z 
Dy (§ (6), (0) 
Z (t,o) 
FB1, z=s (t) or 
t = w (2) 
0 


Fig. 3. Characteristic solution domain for the case g'(t)=0, 
O<t<T 


A 
FB3 FB2 
D2 
P (z ate ) 
D, FB1, z=s (t) or 
t=w (z) 
z = 2 (t,o) (§ (9), (o)) 
fox 
4 


Fig. 4. Characteristic solution domain for the case when 
e=constant > 0, -e<g'(th<0,0<t<T 


It is seen from Eq. (56) that the maximum distance that 
the characteristics will travel is given by 


K, (9,—9%\" 
s u 0 
which is the same as given by Eq. (45). This shows that 


FB2 (the locus 6(z, t) = 9,) is a vertical wall with lower 
end given by Eq. (45) and the upper end by Eq. (56) and 


=T+—(6,— 64). (61) 


In the case when g(t) is not constant, we get the paramet- 
ric representation of FB2: 


K — 0, |* 
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DRAINAGE PERIOD = 6 H, e = 0. 


BROOK AND COREY RELATION 
FIXED DEPTH = 10, 20, 30 CM 


MOISTURE CONTENT 


2/DAY 


Fig. 5. Moisture content as a 
function of time for fixed 


values of depth 


DRAINAGE PERIOD = 6 H, e = 0.2/DAY 
BROOK AND COREY RELATION 


FIXED TIME = 2, 7, 10H 


MOISTURE CONTENT 


Fig. 6. Moisture content as a 
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function of depth for fixed 


T 
50 60 
DEPTH, CM 


tep2 + — 9 
e 
Thus 


dz z' (a) _ 4K, (g(a) — 90)" 
dt 


values of time 


Since g’(c) < 0, we have for FB2, —e < g'(a) < 0 > dz/dt 
<0, and g'(a) < —e dz/dt > 0. Under the condition 
—e<g'(a)<0,0<0<T, the front wall of water recedes 
(dz/dt < 0) when its moisture content is 0). In this case, 
the appearance of FB2 is as shown in Fig. 4. If g’(a) = 0, 
0<o<T, the front wall of water is stationary (dz/dt = 0), 
as shown in Fig. 3, when it has moisture content 4). 
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MOISTURE FLUX 


DRAINAGE PERIOD = 6 H, e = 0.2/DAY 
BROOK AND COREY RELATION 
FIXED DEPTH = 10, 20, 30 CM 


Fig. 7. Moisture flux as a 
function of time for fixed 
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MOISTURE FLUX 


values of depth 


DRAINAGE PERIOD = 6 H, e = 0.2/DAY 
BROOK AND COREY RELATION 
FIXED TIME = 2, 7, 10H 


Fig. 8. Moisture flux as a 
function of depth for fixed 


DEPTH, CM 
Total moisture content 


The total moisture content W varies within the region, 
and can be determined by integrating the moisture profile 
from z = 0 to a particular depth z = L. Let L < z*. Refer- 
ring to Fig. 1, at t = 0, the initial total moisture content 
W,=0,L. For 0<t<T, say, t=t,, when g(t) =con- 
stant, 


W = 0,z,+ 0,(L—z;,) (65) 


values of time 


where z, is the point of intersection of t = t, and z=s(t). 
For t > t*, 


W(t, L) = t) dz (66) 
0 


where 6(z, t) is prescribed by Eq. (55). An explicit solution 
of this integral is, however, not tractable. 
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DRAINAGE PERIOD = 6 H, e = 0.2/DAY 
BROOK AND COREY RELATION 
FIXED DEPTH = 10, 20, 30 CM 
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AVERAGE VELOCITY, CM/H 


Fig. 9. Average velocity as a 
function of time for fixed 
values of depth 


DRAINAGE PERIOD = 6 H, e = 0.2/DAY 
BROOK AND COREY RELATION 
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Fig. 10. Average velocity as 
a function of depth for fixed 
values of time 


DEPTH, CM 


197 
5 
4 
‘| 
il 
it 
tl 
' | 
| 
10 cm 
2 TY 20 cm 
30 cm 
wt 
| \ 
| 
\ 
0 10 20 30 
TIME, H 
5 
4 
{ 
| 
t 
| 
2 10 h 
| H 
q 
| 
/ 
H 
0 10 20 30 40 50 60 70 80 90 100 


198 


An example of application 


A wetting event with moisture applied at the ground sur- 
face was considered, with K (0) given by Eq. (4). The soil 
was Glendale clay loam, investigated by Sisson et al. 
(1980), with the following properties: 0,=0.52, 0) = 
0.246, K,=1m/day, and a= 4.25. The moisture was 
applied at the rate of 2 cm/h for 6 hours. These values 
were used by Charbeneau (1984). The plant-root extrac- 
tion e(t) was considered constant but values of e = 0.1, 
0.2, 0.3, 0.4, and 0.5/day were considered. It should be 
emphasized that these values were used to illustrate the 
procedure presented in the paper. If n=0.52, then 
e =0.1/day would correspond to r(t) = 0.052/day. For 
one meter of soil column, this would translate into 
5.2mm of plant-root extraction for a day. Likewise, 
e = 0.5/day would yield r(t) = 0.25/day or 25 mm/day. 
This rate is clearly high but the objective here was to 
intentionally show the modulation of soil moisture profile 
by the plant-root extraction. 

The kinematic wave model was applied to compute 
A(z, t), q(z, t), and u(z, t). Then @(t), u(t) and q(t) were 
obtained at selected values of z. Likewise, 0(z) and q(z) 
were computed at selected values of t. For a sample case 
of e=0.2/day, computations are illustrated. Figure 5 
shows 0(t) at z = 10, 20, and 30 cm, and Fig. 6 shows 0(z) 
at t=2, 7 and 10h. The moisture profile has realistic 
appearance. The moisture flux is plotted as a function of 
t in Fig. 7, and as a function of z in Fig. 8. The average 
velocity is plotted as a function of t in Fig. 9, and as 
function of z in Fig. 10. Both q and u appear to be realistic 
in their variation. The effect of increasing e was to com- 
press domain D,. In other words, after the moisture flux 
was ended at the ground surface, the moisture profile 
returned to its initial state more quickly where the mois- 
ture extraction rate was higher. For higher e, the charac- 
teristics in domain D, were more curved. 


Conclusions 


The kinematic wave model for soil moisture movement 
has been formulated, and its solution domain involves 
three free boundaries. Analytical solutions are tractable 
for the case when the upstream boundary condition is 
independent of time. For a time-dependent boundary 


condition, numerical solutions are the only resort. In this 
model, the advancing front of water or wetting front is 
necessarily an advancing shock wave. It is plausible to 
assume that the moisture content at the wetting front is 
@,. Since the kinematic wave theory cannot account for 
this assumption, the diffusion wave model has to be em- 
ployed in that case. 
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Abstract. The threshold canopy temperature method for 
controlling a drip irrigation system includes a physiolog- 
ically based threshold temperature and irrigation appli- 
cation rate that responds to the environment. Energy 
input from the environment causes canopy temperature 
to exceed the threshold value and irrigation is then ap- 
plied. This study evaluated temperature distributions, 
amount of optimum time, and the amount of irrigation 
control time for cotton where irrigation scheduling was 
controlled by different threshold temperatures during the 
years 1988 to 1991. Optimum time for cotton growth was 
defined as the accumulated time that canopy tempera- 
tures were between 25 and 31 °C and the time accumulat- 
ed above different threshold temperatures was designated 
as irrigation control time. Threshold temperatures over a 
26 to 32°C range altered the frequency distribution of 
temperature within the optimum temperature range (25— 
31°C) by reducing temperatures above the threshold. 
Frequency of canopy temperatures of a 28°C threshold 
temperature treatment decreased in the 28 to 29°C incre- 
ment and then remained below air temperature. Irriga- 
tion control time was more sensitive than optimum time 
to changes in threshold temperature between 26 and 
31°C. Optimum time and irrigation control time of the 
28°C threshold temperature varied by 37% and 29%, 
respectively. Lint yields in 1988 and 1990 were high while 
those in 1989 and 1991 were low because of unfavorable 
weather. Irrigation amounts applied during DOY 198- 
273 that were above 20 cm in high yield years or 12 cm in 
low yield years did not increase yield. 


Temperature significantly influences plant vegetative 
growth and phenological development. The metabolism 
of plants is affected by temperature in such a manner that 
there is an optimal thermal range for their growth and 
development (Burke et al. 1988). The existence of an op- 
timal thermal range implies that some plant temperatures 
are more conducive to plant growth than others and that 
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an increase in the amount of time that the plant is within 
the optimal range will result in increased plant growth. 

Continuous measurement of canopy temperature and 
application of irrigation for 15-minute periods when a 
threshold temperature was exceeded has been demon- 
strated as an automated irrigation scheduling procedure 
for cotton using a drip irrigation system (Wanjura et al. 
1992). Plant temperature is influenced by its water status 
(Idso et al. 1982; Grimes et al. 1987) which can be altered 
through irrigation. Two assumptions that are made in 
using canopy temperature for irrigation management are 
that, first, canopy temperatures above the threshold indi- 
cate inadequate plant water status and second plant 
growth or productivity are positively correlated with the 
amount of time that the plant’s temperature is within its 
optimal thermal range. 

Optimum time for crop production is a concept that 
recognizes there is a specific temperature where growth is 
maximized. A field grown crop does not grow under a 
constant temperature regime but rather in a fluctuating 
thermal environment. Burke et al. 1988 defined the ther- 
mal kinetic window for cotton as 23.5—32°C. Earlier re- 
search by Arndt (1945) found that the optimum tempera- 
ture for cotton growth was 27°C after germinating for 
three days. An analysis of the thermal dependence of 
reaction velocity of the enzyme glutathione reductase 
from cotton showed that velocity was highest at 28°C. In 
this paper an optimal thermal range of 25—31°C was 
selected where reaction velocity is predicted to be within 
33% of its maximum value. Thus the optimum time 
(OPT) for cotton is the time that the canopy temperature 
is within a range from 25—31 °C, Fig. 1. Mathematically, 
OPT = < time, when 31 °C < T, > 25°C. Optimum time 
can be determined from measured canopy temperatures 
or specified by values of environmental factors that allow 
the plant canopy to maintain its optimal thermal range. 
Since plant canopy temperature depends on both the at- 
mospheric environment and the water status of the plant, 
a well-watered plant was used as a reference. 

In temperature-controlled-irrigation scheduling (TCIS) 
using a drip irrigation system the threshold temperature 
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determines when and how much water is applied. The 
cumulative time that the threshold temperature is exceed- 
ed represents the maximum available time for irrigation 
and is termed irrigation control time (ICT), Fig. 1. Ex- 
pressed mathematically as ICT = 2 time, when T, > 7,, 
where 7, is canopy temperature and 7, is the threshold 
canopy temperature level. Irrigation control time de- 
pends on the threshold canopy temperature level and the 
environment. The overall result is that the threshold tem- 
perature level determines irrigation control time which 
controls the amount of water applied. Applied water af- 
fects the range and distribution of canopy temperature. 

The TCIS method assumes that plant canopy temper- 
ature is a sufficient indicator of plant water stress, see 
Fig. 1. Canopy temperature is a direct indicator of crop 
water status and atmospheric energy dynamics. Canopy 
temperature is used to control irrigation by automatically 
applying irrigation through a drip system whenever a 
threshold canopy temperature is exceeded. The use of 
TCIS requires that the energy input from the environ- 
ment be large enough to raise the crop canopy tempera- 
ture above the threshold temperature. 

In order to assess the effect of canopy temperature 
changes that result from TCIS, an estimate of the maxi- 
mum amount of time that the canopy temperature could 
be within a specific thermal range is required. Drip irri- 
gated cotton where water application was controlled by a 
threshold temperature of 28°C produced the highest av- 
erage yield during a four-year study (Wanjura etal. 
1992). The thermal environment of cotton where irriga- 
tion was controlled by a 28°C threshold is used as a 
reference because yield was maximum and presumably 
water status was optimized. The primary motivation for 
analyzing the thermal environment of the canopy is to 
determine whether the accumulation of temperature 
above a threshold limits the use of threshold canopy tem- 
perature to control irrigation scheduling. 


The purpose of this paper is to examine the tempera- 
ture environment of drip irrigated cotton where schedul- 
ing was controlled by threshold temperature. The specific 
objective was to evaluate temperature distributions, 
amount of time within the optimum temperature range, 
and irrigation control time. 


Materials and methods 


The analysis was conducted using data from a cotton production 
area where irrigation supplements rainfall and the length of the 
growing season is often limited by low temperature. The period of 
the growing season that was used in the analysis was the middle 
portion of the crop production cycle (DOY 198-273) when fruit are 
initiated and matured, and most irrigation is applied. Plant canopies 
were the maximum size during this period which was optimum for 
measuring canopy temperature with infrared thermometers. 

The data was collected from TCIS studies conducted on the 
Texas Agricultural Experiment Station located at Lubbock, Texas 
between 1988-1991. The experimental plots were flat-broken with 
a moldboard plow each year prior to planting and herbicide broad- 
cast and incorporated. Fertilizer application, based on soil samples, 
was applied before planting. A preplant furrow irrigation was ap- 
plied when needed to ensure adequate moisture for seed germina- 
tion and to refill the top 2 m of soil to field capacity. 

Cotton was planted in rows oriented east-west spaced 0.76 m 
(30 inches) apart. Plot size was 18 rows wide by 30.5 m long. Temik 4 
was applied in the seed furrow at planting for the control of thrips. 
Irrigation was applied through 16mm diameter polyethylene 
dripline emitter hose with 61cm spacing between emitters. The 
dripline hose was placed on the soil surface in each row. Nitrogen 
was applied to the plots through the irrigation system at the rate of 
0.221 kg(N)/ha-mm. 

All plots were instrumented as follows: Infrared thermometers 
(IRTs) were placed on each side of a row to view the upper one-third 
of the canopy, and net radiation was measured above the canopy. 
One anemometer and pyranometer were positioned 2 m above the 


' The use of trade names does not imply endorsement by the U.S. 


Department of Agriculture 


Threshold Temperature Controlled Irrigation Scheduling 


Irrigation 


Canopy 
Temperature 


Optimum 
Temperature 
(25 to 31°C) 


®< 
A 


Threshold 
Canopy 


Temperature 


«Optimum 
Time 


Irrigation 
* Control Time 


Fig. 1. Representation of threshold temperature 
controlled irrigation scheduling which depicts the 
relationship of optimum time and irrigation control 
time with specific parameters 
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soil surface in a well-watered treatment to measure ambient condi- 
tions. Dry bulb air temperature (AT) and wet bulb air temperature 
were also measured at a 2 m height in the study area. Instruments 
were located in the north-east quadrant of each plot to provide 
maximum fetch for a prevailing south-west wind direction. 

Irrigation scheduling of threshold canopy temperature con- 
trolled treatments was controlled by two 4° field-of-view IRTs view- 
ing the canopy perpendicularly from a 45° angle from the north and 
south sides. Additional information on temperature measurement 
and data recording are reported in Wanjura et al. (1992). 

The following threshold canopy temperatures were used: 28, 30, 
and 32°C in 1988; 28°C in 1989; 26, 28, and 30°C in 1990; and 28°C 
in 1991. The treatment whose irrigation was scheduled by a 28°C 
threshold temperature (T 28C) in each year was used as the refer- 
ence ‘“‘well-watered” treatment. 

Amount of irrigation control time for different canopy temper- 
ature thresholds were calculated. Optimum time was the amount of 
time that temperatures were between 25—31 °C. Irrigation-control- 
time was the amount of time that canopy temperature exceeded 
different thresholds. 


Results and discussion 
Air and canopy temperatures 


Frequency distributions of dry bulb air temperature at 
2 m (AT) and canopy temperatures of the 28 °C treatment 
(T 28 C) measured during 24-hour periods are compared 
in Fig. 2. The most frequently occurring temperatures in 
all years were 19°C and 18°C, respectively, for AT and 
T28C over the period of day-of-year (DOY) 198-273. 
Above 19°C, AT frequencies declined linearly to an upper 
limit temperature of 33°C. Canopy temperature frequen- 
cies declined from 18°C to 20°C, stabilized between 20 
and 28°C, and declined again between 28°C and 32°C. 
Presumably transpiration cooled the canopy and main- 
tained uniform temperatures in the 20 to 28 °C range. The 
steep decline above 28°C coincides with the threshold 
canopy temperature of 28°C which caused irrigation to 
be applied and reduced the occurrence of temperatures 
above the threshold value. 

The frequency distribution of AT and T 28C was also 
compared in the optimum temperature range between 
25—31°C. The daily period for making the comparison 
was restricted to those times when AT was 25°C or higher 
and total solar radiation was above 200 W/m?. These 
restrictions represent atmospheric conditions when well- 
watered plants were transpiring in response to evapora- 
tion demand. The amount of time that AT and T28C 
were in the optimum range was lowest in 1989, Fig. 3. Air 
temperatures between 25 and 31 °C were more uniformly 
distributed than T28C. The difference (T28C-—AT) in 
daily time that temperatures were within 1 °C increments 
between 25°C and 31°C were statistically compared us- 
ing a 7-test (analysis not shown). These comparisons 
were significant, except when designated as “ns” in 
Fig. 3. 

The effect of using the 28°C threshold canopy tem- 
perature to apply irrigation is shown by comparing the 
cumulative times that T28C and AT were between 
25-29°C and 29-31°C, Table 1. Cumulative time at 
temperatures from 25 through 28°C for T28C was al- 
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Fig. 2. Frequency distribution of air temperatures (AT) and canopy 
temperatures of the 28°C threshold treatment (T 28 C) during 24-h 
periods for the interval DOY 198-273, 1988-1991 


ways greater than for AT; however, from 29 to 31°C AT 
was always greater than T 28C. Irrigating when canopy 
temperature exceeded 28°C reduced the amount of time 
that the cotton canopy exceeded 28°C compared to that 
for AT. Among years both AT and T28C were lowest in 
1989. The only year effect on AT and T 28C was a lower 
amount of time that these temperatures were in the range 
from 25—31°C in 1989. 


Optimum temperature 


Optimum time was calculated using the time that either 
T28C or AT was within the range of 25—31 °C during the 
entire 24-hour period without filtering for any environ- 
mental condition. The distributions of AT and T28C 
within the optimum temperature range were different in 
all years, Fig. 4. The linear relationships between cumu- 
lative frequency with increasing temperature show that 
within the optimum temperature interval AT tempera- 
tures are uniformly distributed. Cumulative frequency of 
T28C increased most rapidly through the 27—28°C in- 


crement but then decreased in comparison to AT. Times 
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Fig, 3. Average daily time that air 
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temperature (AT) and canopy 
temperature of the 28°C threshold 
treatment (T 28 C) were within 
1°C intervals in the optimum 
temperature range (25—31°C) 

for the interval DOY 198-273, 
1988-1991. The daily periods 
were restricted to conditions when 
AT exceeds 25°C and total radia- 
tion exceeds 200 W/m?. Time per 
day is not statistically different 
where “‘ns” appears for the tem- 
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Fig. 4. Comparison of cumulative frequencies of air temperature 
(AT) and canopy temperature of the 28°C threshold treatment 


(T28C) within the optimum temperature range for the conditions 


described in Fig. 3 


at optimum temperature were similar for AT and T28C 


in each year except 1989 during DOY 198-273, Table 2. 


The largest difference in optimum time between AT and 


the T28C occurred in 1989 when AT was 109 h higher 
than T28C. The large difference in 1989 contrasts with 
all other years when differences were 30 h or less. 


Irrigation control time 


Irrigation control time was analyzed by comparing the 
different threshold temperature levels used to schedule 
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Table 1, Cumulative time that canopy temperature of well-watered 
cotton (T 28C) and air temperature measured at two meters (AT) 
were in the temperature range 25—31°C while total radiation ex- 
ceeded 200 W/m? and AT was equal to or greater then 25°C during 
the period DOY 198-273 


Total 
Time* 


Year Days 25-—29°C 29-—31°C 


F28C AT T28C AT 


Hours” 


1988 
1989 60 
1990 62 
1991 55 


66 461 354 
200 14 
384 286 


378 311 


203 38 
121 17 
239 43 
214 37 


123 
49 


90 
106 


* Total number of hours when AT exceeded 25°C and total radia- 
tion was above 200 W/m? 

> Cumulative times for T28C and AT within each temperature 
range were significantly different except in 1989 


irrigation. The year effect on irrigation control time for 
well-watered conditions was estimated by comparing val- 
ues for T28C, Table 3. Irrigation control time had a 
range of 35h (37%) compared with a range of 107h 
(29%) in optimum time. The 37% range in irrigation 
control time would result in the same proportional range 
of application for a drip system with a constant irrigation 
rate. 

The influence of threshold canopy temperature level 
on irrigation control time was estimated by comparing 
the different temperature thresholds used in 1988 and 
1990, Fig. 5. Variations in irrigation control time among 
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Fig. 5. Comparison of total optimum time and irrigation control 
time for threshold canopy temperature levels in 1988 and 1990 for 
the period DOY 198—273 


Table 2. Optimum time (25—31 °C) for canopy temperature of well- 
watered cotton (T28C) and for air temperature measured at two 
meters (AT) for the period DOY 198-273, 1988-1991 
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Fig. 6. Frequency distributions within the optimum temperature 


range (25-31 °C) for different threshold temperature levels and air 
temperature (AT) for the period DOY 198—273, 1988 and 1990 


Table 3. Irrigation contro] time and optimum time for a 28°C 
threshold canopy temperature controlled irrigation treatment 
(T28C) for the period DOY 198—273, 1988-1991 


Year T28C AT T28C AT 


Hours Percent of total 


1988 453 471 26 
1989 366 475 26 
1990 378 408 22 
1991 437 414 23 
Average 409 442 24 


* Percentages were calculated based on 1824 total hours (76 days 
x 24-hours per day) during DOY 198 and 273 


threshold temperature levels were larger than their corre- 
sponding optimum times. In 1988 threshold temperatures 
of 28, 30, and 32°C had a range of 54h in irrigation 
control time which corresponded with a range of 32 h for 
optimum time. In 1990 irrigation control time of 
threshold temperatures 26, 28, 30°C had a range of 196 h 


Year Irrigation (Hours) 


optimum time” 


1988 102 453 
1989 95 366 
1990 109 378 
1991 130 437 


control time* 


4 Time above a threshold canopy temperature of 28°C 
> Time within the temperature range from 25-—31°C 


while optimum time had a range of only 14h. Threshold 
temperatures between 26 and 32°C were inversely related 
to irrigation control time. Higher threshold temperatures 
results in significantly less irrigation control time than 
lower temperatures but their optimum time was not re- 
duced. Since optimum time was similar among different 
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Fig. 7. Relationship of lint yield with optimum time, irrigation con- 
trol time, and amount of irrigation during the period DOY 198 —273 
for high yield years (1988/1990) and low yield years (1989/1991) 


threshold temperature levels but canopy temperatures 
were modified above the threshold temperature used for 
irrigation control it follows that the distribution of 
canopy temperatures within the optimum interval should 
be different among threshold temperatures, Fig. 6. In 
1988 and 1990 threshold temperatures above 28°C had 
few temperatures below 28°C and more temperatures at 
or above 28°C than did the lower threshold tempera- 
tures. The distributions of AT were more uniform than 
for any threshold canopy temperatures. 


Lint yield 


Well-watered irrigation treatments were used to minimize 
water stress effects on yield, Fig. 7. In this comparison, 
well-watered included threshold temperatures of 26°C 
and 28°C and a treatment that fully replaced soil water 
used on a weekly interval as measured by a neutron scat- 
tering technique. An obvious difference in years is appar- 
ent where 1988 and 1990 yields were higher than those in 
1989 and 1991. Optimum time and irrigation control time 
had similar ranges in both the high and low yield years 


which suggests that they were not the causes for these 
large yield differences. The year effect on yield occurred 
outside the period, DOY 198-273, for which optimum 
time and irrigation control time were accumulated. The 
yield reduction in 1989 and 1991 resulted from plant 
damage caused by cool temperatures or thunder storms, 
(prior to DOY 198) that reduced plant vigor and popula- 
tion and late season cool temperatures that delayed ma- 
turity. The linear and second-order curves were comput- 
ed using regression analysis. Criterion for choosing the 
best fit curve was that all regression model coefficients be 
statistically significant. Optimum time was linearly relat- 
ed to lint yield in high yield years (1988 and 1990) and 
nonlinearly in low yield years (1989 and 1991). Optimum 
time had the greatest effect on yield in the low yield years 
because low temperatures in 1989 and 1991 after DOY 
243 did not allow the crop to fully mature. Irrigation 
control time was linearly related to yield in the low yield 
years and nonlinearly in the high yield years. Again, the 
effect of irrigation control time was greater in the low 
yield years because low temperatures after DOY 243 did 
not permit full crop maturity. 

Total irrigation applied to the well-watered treat- 
ments, in direct proportion to differences in irrigation 
control time, during the period DOY 198-273, had a 
neutral effect on lint yield in the high yield years and a 
negative effect in the low yield years, Fig. 7. In the high 
yield years well-watered treatments received adequate ir- 
rigation during DOY 198-273 with some treatments be- 
ing excessively irrigated. In low yield years irrigation 
amounts that showed a neutral effect in the high yield 
years decreased yield. Since one of the differences be- 
tween high and low yield years was a shortened growing 
season, it is a reasonable assumption that the higher irri- 
gation amounts applied in low yield years slowed crop 
maturity and reduced yield. 


Conclusions 


Correct application of threshold temperature controlled 
irrigation scheduling requires the crop to be thermally 
adapted to its growing environment and the threshold 
temperature must be appropriately selected to result in an 
amount of irrigation (water stress level) that is compati- 
ble with efficient production. The 28 °C threshold canopy 
temperature which results in a well-watered crop condi- 
tion had maximum daily irrigations that ranged from 
10-12 mm during the four years of this study. These 
irrigations were applied through a drip irrigation system 
where the average application rate varied from 1.5 to 
2.1 mm/h. Daily potential evapotranspiration estimated 
by the Penman-Monteith combination method for a well 
watered short grass surface described in Jensen et al. 1990 
averaged 7 mm/day during the period DOY 198-273. 
Thus, the irrigation system and scheduling procedure had 
the capability of completely supplying the full evapotran- 
spiration demand. The irrigation scheduling method used 
in this study resulted in frequent irrigations. Small fre- 
quently applied irrigations are an effective and efficient 
method of irrigation (Bordovsky et al. 1992; Radin et al. 
1989). 
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Different threshold canopy temperature levels did not 
significantly change the amount of optimum time when 
canopy temperature was within the optimum tempera- 
ture range from 25—31 °C, but the frequency distribution 
of canopy temperatures within the optimum range was 
altered. Controlling irrigation scheduling with different 
threshold temperature levels changed the amount of irri- 
gation control time and quantity of irrigation. Irrigation 
control time and optimum time were positively related to 
cotton lint yield. 


References 


Arndt CH (1945) Temperature-growth relations of the roots and 
hypocotyls of cotton seedlings. Plant Physiol 20: 200-220 
Burke JJ, Mahan JR, Hatfield JL (1988) Crop specific thermal 
kinetic windows in relation to wheat and cotton biomass pro- 

duction. Agron J 80:553—556 


205 


Bordovsky JP, Lyle WM, Lascano RJ, Upchurch DR (1992) Cotton 
irrigation management with LEPA systems. Trans ASAE 
35:879-884 

Grimes DW, Yamada H, Hughs SW (1987) Climate-normalized 
cotton leaf water potentials for irrigation scheduling. Agric Wa- 
ter Mgt 12:293-—304 

Idso SB, Reginato RJ, Reicosky DC, Hatfield JL (1982) Determin- 
ing soil-induced plant water potential depressions in Alfalfa by 
means of infrared thermometry. Agron J 73:826—830 

Jensen ME, Burman RD, Allen RG (1990) Evapotranspiration and 
irrigation water requirements, ASCE Manuals and Reports on 
Engineering Practice No. 70, American Society of Civil Engi- 
neers, New York, NY, 332 pp 

Radin JW, Mauney JR, Kerridge PC (1989) Water uptake by cotton 
roots during fruit filling in relation to irrigation frequency. Crop 
Sci 29: 1000-1005 

Wanjura DF, Upchurch DR, Mahan JR (1992) Automated irriga- 


tion based on threshold canopy temperature. Trans ASAE 
35:153-159 


| 
| 

j 


q 
i 
q 
q 
| 
q 
4 
a 
3 
| 
| 
| 
| 
| 
| 
| 
i | 
q 
| 
| 
4 


Irrig Sci (1994) 14:207—212 


Irrigation 


Science 
© Springer-Verlag 1994 


Control of irrigation amounts using velocity and position 
of wetting front 


B. Zur', U. Ben-Hanan’, A. Rimmer ', A. Yardeni * 


' Agricultural Engineering Department, Technion, Haifa, Israel, 32000 
2 Mechanical Engineering Section, Ort-Braude College, Carmiel, Israel 


3 Agroteam Consultants Ltd., Migdal HaEmeck, Israel 


Received: 29 May 1993 


Abstract. A new approach for the estimation and control 
of the quantity of water applied in an irrigation is pre- 
sented in which irrigation is stopped when the wetting 
front reaches a critical depth, Z,. An expression for cal- 
culating the critical depth Z, was developed. A major 
parameter in this expression is the velocity of advance of 
the wetting front, V, which was shown to be directly relat- 
ed to the application rate, IR, and inversely related to the 
initial soil water content, 0;. A depth probe (patent pend- 
ing) was designed, constructed and tested for the purpose 
of monitoring the position of the wetting front during 
infiltration and redistribution and for computing the val- 
ue of V. Equations developed for relating the velocity of 
advance of the wetting front to 0; as well as for estimating 
the value of the critical depth Z, were successfully tested 
under conditions of uniform distribution of the initial soil 
water content. An iterative learning process which uti- 
lizes the real time output from the depth probe during 
each irrigation and is therefore capable of handling real- 
istic field conditions where nonuniformity is the rule is 
presented. The acquired information is used to estimate 
a critical depth of the wetting front, Z,, for a planned 
final wetted depth, Z,, during each irrigation. This pro- 
cess is incorporated in the depth probe and is used to stop 
irrigation and thus control the quantity of water applied. 


Ideally the net amount of water applied during an irriga- 
tion, Q, is defined as: 


(1) 


where 0, is the volumetric soil water content at field ca- 
pacity, 0; is the volumetric soil water content prior to 
irrigation and Z, is the desired depth of soil to be wetted 
to field capacity. The value of Z, is estimated using the 
current rooting depth of the crop as well as the leaching 
requirements. 
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An accurate estimate of Q, is a basic requirement for 
optimal irrigation management. Applying the correct 
amount of water is important both for efficient water use 
as well as to reduce groundwater contamination resulting 
from deep percolation. The use of Eq. (1) requires a de- 
pendable estimate of 6, at a range of soil depths, whereas 
the values of 0, at these depths should be known as part 
of the general information representing the soil. Taking 
soil samples to estimate 6; is seldom practiced because it 
is cumbersome and labor consuming. Instead, physical 
and empirical relationships employing meteorological 
parameters and crop coefficients are generally employed 
to estimate water use by the crops. Current procedures for 
estimating Q, are incapable of identifying the depths from 
which soil water was extracted by the crop and lack the 
means to insure that the estimated Q, is sufficient to wet 
the desired soil depth. 

In the present paper we propose a new way to estimate 
the value of Q, and to monitor is distribution in the soil 
profile. In principle, the value of Q, applied in the field is 
determined by the duration of the irrigation. We propose 
to stop an irrigation, and thus control Q,, when the wet- 
ting front reaches a critical depth, Z,. The critical depth 
Z, is defined as the position of the wetting front at the end 
of an irrigation from which sufficient water would be 
redistributed to bring the soil profile down to the depth 
Z, to its field capacity. 

The specific objectives of the present paper are: 

1. Develop equations for estimating the critical depth Z,. 
2. Present a system for tracking the position of the wet- 
ting front during an irrigation. 

3. Test objectives 1 and 2. 

4. Discuss irrigation control algorithms under realistic 
field conditions. 


Theory 


The process of wetting a soil profile goes through two 
stages. 
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At the end of the first, infiltration stage, a soil depth, 
Z,, is wetted in excess to a water content, 6,, which is a 
function of the application rate and is constant with 
depth, Fig. 1. During the second, redistribution stage, ex- 
cess water from the initially wetted soil depth drains 
downward to wet an additional depth of soil. The result- 
ing wetted soil profile down to the depth Z, is character- 
ized by a relatively stable soil water content of field capac- 
ity, 0;, Fig. 1. Duration of the redistribution stage varies 
depending on soil texture, soil structure and depth of 
wetted profile. For medium textured soils the duration of 
this stage is about 24 h. 

Ideally, the arrival of the irrigation water to the depth 
Z,, could be used as a signal to stop an irrigation. Howev- 
er, since the irrigation water reaches the depth Z, after a 
period of redistribution this possibility is not practical. 
An alternative approach is to define a critical depth Z, as 
the position of the wetting front at the termination of 
irrigation which would bring the soil profile down to the 
planned wetting depth Z, to its water content of field 
capacity at the end ofa redistribution period. The arrival 
of wetting front to the depth Z, could be used as a signal 
to stop the irrigation. The relationship between the 
planned depth, Z,, and the critical depth, Z,, is expected 
to depend on soil type, initial soil water content distribu- 
tion and the application rate, ] R. 

Values of Z, which would yield desired depths Z,. fol- 
lowing redistribution could be computed using a numeri- 
cal procedure where the Richards equation is solved for 
the infiltration and redistribution processes (Van Ge- 
nuchten, 1978, 1980). However, this procedure requires 
initial conditions which include knowledge of the initial 
distribution of volumetric soil water content, 0;, which is 


Volumetric Soil Water Content - © 
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Fig. 1. Schematic presentation of water balance during irrigation 
and redistribution 


usually not available. In addition, hysteresis plays a sig- 
nificant role in water flow during redistribution and nu- 
merical solutions are not properly equipped to take this 
phenomenon into consideration (Van Genuchten, 1978, 
1980). Therefore, this possibility is of limited value when 
applied to field conditions. 

In the present paper a simple relationship between Z, 
and Z, is developed based on the equality of the mass 
balance of soil water at the end of an irrigation and after 
24 h of redistribution. For the conservation of mass the 
areas ABDC and EBGF, Fig. 1, must be equal. The area 
ABDC which represents the amount of water added to 
the soil during an irrigation is estimated as: 


Q:= (4, —9;) (2) 


The arrival of the wetting front to the depth Z, is indicat- 
ed in the proposed approach by the response of a sensor 
positioned at that depth. The value of 6, is estimated by 
taking soil samples down the wetted profile immediately 
following an irrigation. 

The area EBGF represents the distribution of the add- 
ed water at the end of redistribution, therefore: 


QO, = (0,—9;) Ze (3) 


The arrival of the redistribution front to the depth Z, is 
indicated in the proposed system by a sensor positioned 
at that depth which had just sensed the arrival of the 
redistribution front about 24 to 48 hours after an irriga- 
tion. Equating Eq. (2) and (3) we obtain: 


(0,—9;) Z, = Ze. (4) 
Similarly, from Fig. 1, the areas AEHC and HDGF are 
equal and we can write: 

(0, Z, AZ (5) 


where the left side of Eq. (5) represents the amount of 
water drained from the initially wetted profile and the 
right side of Eq. (5) represents the additional depth AZ, 
wetted to field capacity. 


Based on Eq. (4) and (5) we can now write: 
(0, Z, (6,—0;)* Zp 
For ideal systems where 6;, 0, and 0, are uniform for 
depths larger than Z,, Eq. (6) becomes: 
AZ (0, — 
7 
Ze 


and since AZ = Z,—Z,: 


Zi (0, — 9) 
®) 


The expression (6, — 0;) represents the volume fraction of 
water in the wetted soil depth which will drain to depths 
below Z, during redistribution. Generally this quantity is 
characteristic of a soil, depends on soil texture and tends 
to increase with sand content. The expression on the right 
side of Eqs. (7) and (8) is the ratio between the volume 
fraction of water drained from the wetted profile during 
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redistribution to the volume fraction of water taken up by 
the same depth during irrigation. 

Equations (4—8) are valid for conditions where 6, is 
uniform with depth, for infiltration times longer than a 
critical value after which 6, is constant with depth, and 
for a piston type advance of the wetting front. An expres- 
sion for the velocity of advance of the wetting front, V, for 
the above conditions was proposed by Phillip (1957) and 
later modified by Rubin and Steinhardt (1963) as follows: 


_ IR-K()) 
°) 


where | R, the application rate, is smaller than the stable 
infiltration rate of the soil. 

The value of K (@;) in Eq. (9) is normally much less 
than J R and therefore can be neglected. Values of J R and 
0, are interrelated and are constant for a given irrigation 
system-soil combination. 

Thus, the velocity of advance of the wetting front, V, is 
expected to be inversely related to the value of (6, —9;) 
and therefore to the value of 0;. Velocity of advance of the 
wetting front, V, is a major parameter in the present 
scheme for estimating the critical depth Z, and for con- 
trolling the value of Q,. 

Introducing Eq. (9) into Eq. (8) we obtain: 


IR 


a= Vz 
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By definition, the velocity of advance of the wetting front 
is the change with time in the position of the wetting front 
during irrigation. Thus, the measurement of V must in- 
volve a system capable of monitoring the position of the 
wetting front during irrigation. The same system could 
signal the arrival of the wetting front to the critical depth 
Z,, for stopping the irrigation, as well as for indicating the 
arrival of the redistribution front to the depth Z,, to de- 
termine whether Q, was sufficient. 


Materials and methods 
Description of system 


Velocity of advance of the wetting front, V, is estimated by noting 
the times at which the wetting front reaches predetermined soil 
depths during an irrigation. This is achieved by using a wetting 
depth probe (patent pending), Fig. 2. Sensors which respond to the 
arrival of the wetting front are imbedded at 5 cm intervals along the 
probe. Each sensor consists of a hollow cylinder of a hydrophilic 
porous material with thin annular stainless steel electrodes tightly 
fitted at the lower and upper ends. The spacing between adjacent 
sensors is taken up by a hollow plastic cylinder of the same outside 
dimensions. The upper annular electrode is connected to a central 
metal rod which acts both as an electrical common and as a pole 
which holds all the parts together. The lower annular electrode is 
connected to a separate wire led through the hollow inside of the 
probe to an external device. The electrical resistance between the 
two annular electrodes changes immediately following a change in 
the amount of water absorbed by the hydrophilic porous cylinder. 


Fig. 2. Schematic presentation of the wetting 
ASSEMBLED 


depth probe 
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The electrical resistance or current flow across each sensor along the 
wetting depth probe is measured at frequent time intervals during 
irrigation and redistribution. A microprocessor circuit positioned at 
the top of the wetting depth is programmed to scan the current flow 
across each sensor at predetermined time intervals, record the time 
at which a predetermined change in the meter reading of each sensor 
takes place and log and store the data. 


Experimental studies 


In 1986, a preliminary study using metal oil drums 56 cm in diame- 
ter and 86 cm high were packed with a Hamra loam soil was con- 
ducted. Prior to packing the soil was moistened to initial water 
contents of 10, 12 and 15% by weight. In addition, the air dry soil 
with a water content of 5% was also used. Four drums packed with 
soils at the above 0; values were irrigated by minisprinklers at JR 
values in the range of 7-8 mm/h. An additional group of tests with 
similar 0, values were subjected to 7R values in the range of 26— 
32 mm/h. The velocity of advance of the wetting front as well as the 
values of 0,, 0,, Z,, and Z, were measured using an older version 
of the wetting depth probe. 

In later experiments, plastic drums with a diameter of 54 cm and 
height of 96 cm were packed with soil of known initial water content 
0;, to a predetermined dry bulk density. Two or three wetting depth 
probes were placed in each drum during packing. Precautions were 
taken to assure good contact between soil and sensors during pack- 
ing of the drums. Prior to irrigation soil samples were taken at 
10 cm depth intervals to estimate the value of 6; in the drum. The 
sampling holes were than refilled with the same soil. Mini-sprin- 
klers, with application rates calibrated on site, were used to apply 
controlled amounts of water at known application rates to the soil 
in the drums. Scanning of the sensors on each probe at five minutes 
intervals started 20 minutes before an irrigation and continued at 
this rate during water application. Time intervals for scanning were 
reduced to once every hour when irrigation was terminated. Soil 
samples for water content determination were taken at 10 cm depth 
intervals immediately following an irrigation to estimate the value 
of 0,, and after 24h to estimate 06,. 

The time the wetting front reached a sensor at a given depth was 
obtained from reviewing printouts of the scanning results and not- 
ing the times when the current flow across each sensor changed by 
more than 5% of the original reading. Velocity of advance of the 
wetting front, V, was obtained from linear regression of the relation- 
ship between the time of arrival of the wetting front and depth. 
Values of Z, and Z, were estimated from the same printouts during 
irrigation and redistribution and were compared with estimates 
from the results of the gravimetric soil sampling. 

Three soils were used in the present trials, a sand, a sandy loam 
and a kaolinitic clay. Table 1 summarizes the dry bulk densities, 
initial soil water contents, 0;, application rates, 7/R, maximal soil 
water content at the end of an irrigation, 0,, and soil water content 
after 24h of redistribution, 0,, for the three test soils. 


Results and discussion 


Results of scanning the depth probe every five minutes 
during infiltration into a Hamra sandy loam soil at an 
initial volumetric water content of 3% are plotted in 
Fig. 3. Based on the time when a detectable change in the 
meter readings was observed, it is clear that the wetting 
front passed by the sensor positioned at the 5cm depth 
115 min after the start of the irrigation. Similarly, the 
wetting front passed by the sensor at 10 cm depth after 
195 min, by the sensor at 15cm after 280 min, by the 
sensor at 20 cm after 370 min and by the sensor at 25 cm 
after 460 min. Irrigation was terminated after 505 min. 
The response of the sensors to the arrival of the 
wetting front was instantaneous and the changes in the 


Current ma. 


Time min. 


Fig. 3. Response of individual sensors during infiltration into Ham- 
ra sandy loam soil. Number in box indicates depth in cm 


Table 1. Basic parameters of test soils 


IR 
mm/h 


Soil type 04 0; 
gm/cm* 


Sand 15 0.5 10.0 
6 10.0 
7 10.5 


Sandy loam 3 
4 7.0 
10 9.1 

14 10 


Bukata clay ! 15 10.0 
(kaolinitic) 34 0.77 
36 0.95 47 


* Water drained to bottom of drum during redistribution 


meter readings were considerable. The meter readings for 
each sensor reached a new stable value 15 to 25 min after 
the initial change. The average velocity of advance of the 
wetting front, V, in this trial was computed to be 3.0 cm/h 
with a R? value of 0.98. 

The time required for the wetting front to reach speci- 
fied depths during infiltration into a Bukata Kaolinitic 
clay soil with initial volumetric soil water contents of 15 
and 36% are presented in Fig. 4. Three depth probes were 
used in the drum with 15% volumetric water content and 
only two probes in the drum with soil at 36% water 
content. Computed average velocities of advance of the 
wetting front from the reading of the three probes were 
2.6, 2.7 and 2.7 cm/h for the 15% initial volumetric soil 
water content. The two probes placed in the drum with 
soil at volumetric water content of 36% resulted in aver- 
age velocities of 7.6 and 7.7 cm/h. 

Similar results for a Hamra sandy loam soil with initial 
volumetric water contents of 3 and 10% are presented in 
Fig. 5. For the soil with an initial volumetric water con- 
tent of 3% velocities of 4.5, 4.2 and 4.6 cm/h were comput- 
ed using linear regression. For an initial volumetric water 
content of 10% velocities of advance of the wetting front 
of 10.8, 10.8 and 9.8 cm/h were computed. 
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Fig. 4. Time of arrival of the wetting front as a function of depth 
under two initial volumetric water contents during infiltration. 
Bukata clay soil 
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Fig. 5. Time of arrival of the wetting front as a function of depth 
under two initial volumetric water contents during infiltration. 
Hamra sandy loam soil 


Results of tests with a sand at initial volumetric water 
contents of 0.5 and 6% are presented in Fig. 6. Velocities 
of advance of the wetting front of 4.7 and 5.3 cm/h were 
computed for the sand soil at an initial volumetric water 
content of 0.5%. Under an initial volumetric water con- 
tent of 6% velocities of 6.5 and 6.6 cm/h were computed. 

Results presented in Table 1, in Figs. 4, 5, 6 and from 
additional measurements are used to test the applicability 
of Eq. (9) for estimating the value of 6, from measured 
values of V, IR and 6,. Computed values of 6;, using 
Eq. (9) are plotted against measured values in Fig. 7. The 
points in this figure represent a range of soils, application 
rates and initial soil water contents. A linear regression 
analysis of the data yielded an R? value of 0.98 and a slope 
of 0.98. 
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Fig. 6. Time of arrival of the wetting front as a function of depth 
under two initial volumetric water contents during infiltration. 
Sandy soil 
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Fig. 7. Relationship between measured and computed values of the 
initial soil water content 


The assumption of a piston type advance of the wet- 
ting front can be evaluated from the response of individu- 
al sensors to the passage of the wetting front during infil- 
tration. The response of sensors to a change in the water 
content of the surrounding soil was instantaneous. Thus 
the thickness of the wetting zone can be estimated from 
the product of the time interval between the first response 
of a sensor and when its reading became stable, and the 
average value of V. Based on such an analysis and under 
the present test conditions the thickness of the wetting 
zone for the Bukata kaolinitic clay soil ranged from 2 cm 
when 0; was that of an air dry soil to 4cm when 0; was 
36%. For the sand and Hamra sandy loam soils the thick- 
ness of the wetting zone was in the range of 2—4 cm under 
all test conditions. 
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Fig. 8. Relationship between observed and predicted ratio of rela- 
tive depth of wetting, Z,/Z, 


Table 2. Coefficients from linear regression analysis 


Intercept R? 


1991 data 
1986 data 


data 


—0.063 
0.055 


0.025 


The relationship between measured and predicted Z,/ 
Z,, ratios from both the plastic container and metal oil 
drum experiments representing a range of soil types, ini- 
tial water contents and application rates are presented in 
Fig, 8. Estimated values of the ratio Z,/Z, were comput- 
ed Eq. (10). The ability of Eq. (10) to predict the ratio 
Z,/Zy Was evaluated by subjecting the data presented in 
Fig. 8 to linear regression analysis. Three linear regres- 
sions tests were performed. In the first test only the 1991 
data was used, in the second test only the 1986 results 
were analyzed and in the third test all the data was pooled 
together. The resulting coefficients of the linear regression 
analysis are presented in Table 2. In all three linear regres- 
sion tests the slope of the regression lines was between 3 
to 8% higher than 1.0. The R? value of 0.87 for the 1991 
data suggests a good agreement between computed ratio 
of Z,/Z, using Eq. (10) with the improved wetting depth 
probe and measured values. An R? value of 0.74 was 
obtained for the 1986 data and 0.644 for the pooled data. 
In all three tests the computed intercepts were close to 
zero. Since the sensors were spaced at 5 cm intervals the 
uncertainty in determining the position of the wetting 
front was +2.5 cm. 

When the initial soil water content distribution is uni- 
form down to depths below Z,, the use of Eq. (10) to 
estimate the critical depth of the wetting front Z, was 
shown to be successful. However, under field conditions 
the distribution of 6; prior to an irrigation is far from 
being uniform. Under these conditions the distribution of 
0; below the critical depth, Z,, can not be estimated from 
measurements of the velocity of advance of the wetting 
front during an irrigation. Without this information the 
resulting final wetted depth, Z,, cannot be computed. 


Thus, under realistic field conditions Eq. (10) cannot 
be used to estimate the critical depth Z,. Instead, an 
iterative learning process for estimating the value of Z, 
was developed. The learning process is built on the fol- 
\owing statements which were based on our laboratory 
and field experience. The velocity of advance of the 
wetting front during irrigation is essentially constant and 
is responsive to an average value of the initial soil water 
content distribution and to soil texture. Following termi- 
nation of irrigation the velocity of advance of the wetting 
front gradually slows down and approaches zero at the 
final wetting depth, Z,. The redistribution depth (Z,—Z,), 
is related to the velocity of the advance of the wetting 
front during irrigation. 

The iterative learning process starts during the first 
irrigation of the season. A desired final wetting depth, Z,, 
is introduced and a critical depth Z, is estimated by mul- 
tiplying the value of Z, by a reasonable coefficient in the 
range of 0.4 to 0.7 depending on soil type. Irrigation is 
then initiated and the time of arrival of the wetting front 
to each of the sensors spaced along the depth probe is 
recorded and the average velocity of advance of the 
wetting front computed. Irrigation is stopped when the 
wetting front reaches the sensor closest to the estimated 
depth Z,. The time of arrival of the wetting front during 
redistribution to the deeper sensors and the depth of the 
wetting front 24 and 48 h after stopping the irrigation are 
recorded. Coefficients are than computed and introduced 
into relationships based on the above statements. At the 
start of the next irrigation the probe is activated and the 
average velocity of advance of the wetting front is com- 
puted based on measurements in real time. A new value 
for the critical depth Z, is than computed based on the 
present velocity of advance of the wetting front, on the 
site calibrated empirical relationships and on the planned 
value of Z,.. Irrigation is stopped when the wetting front 
reaches the depth Z,. The position of the wetting front 
during redistribution and the depth of the wetted profile 
24 and 48 h after irrigation are recorded. In case there is 
a significant difference between the planned and mea- 
sured values of Z,, a new value of Z, is calculated and 
tested during the next irrigation. This iterative learning 
process is repeated during each irrigation and changes are 
made in the value of Z, whenever necessary. The de- 
scribed system does not require initial calibration of any 
kind, the sensors respond to relative change rather than 
to quantities or potentials and the decision when to stop 
an irrigation is based on real time measurements. The 
manufactured system is presently under going field tests. 
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